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ABSTRACT 


Comparisons  of  experimentally-extracted  pressure/rate-of-strain  values  were  made 
to  theoretical  pressure-strain  models  for  several  locations  of  a  wing-body  junction  flow. 
The  experimental  pressure/rate-of-strain  results  were  calculated  from  data  obtained  with 
a  LDV  technique.  The  data  consist  of  profiles  of  mean  velocity  and  higher  order  moments 
including  quadruple  products  in  a  two-dimensional  turbulent  boundary  layer  (2DTBL),  a 
strongly  skewed  three-dimensional  turbulent  boundary  layer  (3DTBL),  in  the  vicinity  of  a 
3-D  separation  line,  and  around  the  center  of  the  horse-shoe  vortex  that  forms  around  the 
wing.  Terms  in  the  transport  equations  for  the  Reynolds’  stresses  and  the  turbulent  kinetic 
energy  are  also  presented  here.  Several  linear  and  non-linear  pressure/rate-of-strain 
models  are  tested  using  the  measured  quantities  as  input  to  the  models.  The  tested 
models  are  the  Launder-Reece-Rodi  (1975),  Gibson-Launder  (1978),  Oberlack-Peters 
(1993),  Fu-Launder-Tselepidakis  (1987)  (2  models),  Shih-Lumley/Choi-Lumley  (1985, 
1984),  and  Speziale-Sarkar-Gatski  (1991).  The  measurements  were  carried  out  in  the 
Virginia  Tech  Boundary  Layer  Tunnel,  at  nominal  air  speed  of  27.5  m/s  around  a  NACA 
0020  tailed  and  3:2  elliptical  nosed  wing  shape. 

Near-wall  correction  effects  on  the  theoretical  pressure/rate-of-strain  models;  effect 
of  approximation  for  the  pressure  diffusion,  the  anisotropic  dissipation  approximation  as 
well  as  isotropic  dissipation  approximation  on  the  experimental  pressure-strain  data  are 
examined.  The  results  are  compared  to  previous  Direct  Numerical  Simulation  data 
qualitatively.  Data  show  that  the  pressure-diffusion  approximation  by  Lumiey  (1978)  is  an 

important  term  in  the  ,  W  ,  and  viv  stress  budgets.  Qualitative  comparison  of 

experimental  budgets  of  transport  of  stresses  to  DNS  solutions  show  a  better  agreement 
between  the  two  with  the  anisotropic  dissipation  of  Hallback,  Groth,  and  Johansson 
(1 990).  Pressure-strain  model  estimates  deteriorate  with  increased  complexity  of  the  flow. 
The  models  agree  with  the  data  better  relatively  away  from  the  wall.  Superiority  of  the 
models  depends  on  the  approximations  used. 
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XIV 


CHAPTER  I  -  INTRODUCTION  AND  EXPERIMENTAL  DATA 


The  data  presented  here  are  part  of  a  set  of  measurements  made  in  the  horse¬ 
shoe  vortex  that  forms  around  a  wing-body  junction  vortex.  The  wing  used  is  a  3:2 
elliptical  nosed,  NACA  0020  tailed  wing  profile  with  maximum  thickness  of  2.824  inches, 
chord  iength  of  12  inches,  and  a  height  of  9  inches.  The  wing  was  mounted  perpendicular 
on  the  flat  bottom  surface  of  the  Boundary  Layer  Tunnel  with  a  3  feet  x  10  inches 
rectangular  cross  section  and  8  m  long  tunnel  test  section  at  Virginia  Tech.  The  nominal 
speed  of  air  was  27.5  m/sec  and  Reg  (Reynolds  number  based  on  momentum  thickness) 
of  the  flow  at  0.75  chord  upstream  of  the  wing  was  5936. 

Previous  measurements  using  the  same  tunnel  and  wing  configuration  were 
reported  by  several  researchers.  The  approach  flow  in  front  of  the  wing  at  the  symmetry 
plane  were  reported  by  Devenport  and  Simpson  (1990),  and  Shinpaugh  and  Simpson 
(1995).  Three-dimensional  turbulent  boundary  layer  (3DTBL)  measurements  on  the  tunnel 
floor  due  to  the  presence  of  the  wing  away  from  the  horse-shoe  vortex  were  reported  by 
OIgmen  and  Simpson  (1995a).  Measurements  of  3DTBL  at  the  side  and  at  the  wake  of 
the  wing  were  reported  by  Fleming  et  al.  (1993).  Ha  and  Simpson  (1993)  analyzed  the 
flow  structure  using  several  hot-wire  probes  and  surface  hot-film  probes  at  the  same 
locations  of  OIgmen  and  Simpson  (1995a).  OIgmen  and  Simpson  (1992)  studied  the 
surface  pressure-fluctuations  on  a  stagnation  line  on  the  symmetry  plane  of  the  wing. 
Lewis  et  al.  (1994)  studied  the  heat  flux  fluctuations  at  the  same  locations.  Using  the 
similar  shape  but  a  water  tunnel  many  researchers  studied  the  flow  using  flow 
visualization  methods  (Fleming  et  al.  (1995),  Kim  et  al.  (1991)) 

The  flow  is  well  documented.  The  mean  surface-pressure  distribution  was 
measured  and  reported  by  OIgmen  (1990).  Oil  flow  visualization  methods  were  used  to 
document  the  surface  skin  friction  lines.  In  this  Chapter  several  profiles  measured  by  the 
LDV  technique  using  the  three-component  portion  of  the  five-component-LDV  probe 
(OIgmen  and  Simpson,  1995b)  developed  at  Virginia  Tech  to  study  the  structure  of  the 
3DTBLS  are  presented. 

Figure  1  shows  the  schematic  of  the  wing-body  junction.  The  measurement 
locations  on  the  line  passing  through  Station  5  and  perpendicular  to  the  wing  (23°  line) 
and  the  coordinate  system  definitions  are  also  shown.  Opg  stands  for  the  free-stream 
coordinates  and  the  ()jc  stands  for  the  tunnel  coordinates,  y  axis  is  perpendicular  to  the 
floor.  Positive  x^c  is  along  the  tunnel  axis  looking  downstream  with  its  origin  at  the  wing- 
floor  intersection.  Positive  x  ps  is  in  the  direction  of  the  mean  velocity  direction  at  the  edge 
of  the  boundary  layer.  The  z  axes  complete  right-handed  coordinate  systems.  Figure  2 
shows  the  oil  flow  visualization  picture  obtained  on  the  wall.  Figure  shows  that  incoming 
flow  separates  at  the  nose  region  and  rolls  around  the  wing  to  form  the  horse-shoe  vortex 
region.  At  downstream  of  the  separation  line,  line-of-low-shear  forms.  Line-of-low-shear 
and  the  separation  line  merge  together  downstream  of  the  nose.  Fish-tailed  structure  of 
the  horse-shoe  vortex  is  also  shown.  The  profiles  presented  are  4  out  of  38  profiles 
measured.  These  profiles  presented  were  obtained  in  a  2-D  boundary  layer  without  the 
wing  presence,  and  in  the  horse-shoe  vortex  forms  around  the  wing.  30  profiles  measured 
inside  the  horse-shoe  vortex  were  closely  spaced  in  order  to  calculate  the  x  and  z 
derivatives  of  the  flow  quantities.  Measurements  were  made  on  a  line  passing  through 
Station  5  (OIgmen  and  Simpson,  1995a)  and  perpendicular  to  the  wing  making  =23°  to 


the  z  axis  and  on  a  line  0.125  inches  downstream  of  this  line  and  parallel  to  it.  5  profiles 
were  also  taken  at  Station  5  location.  Measurements  include  the  mean  velocities,  and 
fluctuating  velocity  products  including  4th  order  products.  The  measurements  were  used 
to  calculate  the  terms  in  transport  of  Reynolds  stresses  and  the  extracted  data  were  used 
to  compare  pressure-strain  models  to  the  pressure-strain  data.  Figure  3  shows  the  mean 
velocity  magnitude  and  the  secondary  flow  streamlines  in  the  plane  perpendicular  to  the 
floor  passing  through  the  23°  line.  Abscissa  shows  the  distance  from  the  wing,  where  the 
wing  is  located  at  AS=0.  The  ordinate  is  logarithmically  spaced.  Figure  3  shows  that  a 
large  vortex  is  located  at  the  outer  region  and  entrains  the  high  speed  free-stream  flow 
closer  to  the  wall  especially  closer  to  the  wing.  At  AS=-4.2cm  the  flow  experiences  a 
three-dimensional  separation  where  all  the  secondary  streamlines  merge  towards  a  point. 
At  distances  closer  to  the  wing  (AS»>-0.4cm)  the  flow  direction  is  reversed  (figure  3). 

Measurements  presented  in  this  report  are  mean  velocity,  and  fluctuating  velocity 
products  including  the  3rd  order  at  selected  locations.  Presentation  of  the  other  profiles 
will  be  in  an  accompanying  report.  The  locations  selected  are  the  2-D  station.  Station  5, 
Separation,  and  Vortex-core  stations.  Plots  of  the  measured  profiles  in  a  plane 
perpendicular  to  the  tunnel  floor  encompassing  the  measurement  station  locations  show 
that  the  Separation  and  Vortex-core  stations  correspond  to  the  location  where  the  3DTBL 
goes  through  a  3-D  separation  and  to  the  location  where  the  center  of  the  outer  region 
vortex  is  located,  respectively.  Presentation  of  the  data  and  the  following  discussions  are 
presented  in  tunnel  coordinates  throughout  the  report. 

I-A)  Mean  velocities 

The  U/  mean  velocity  profile  at  the  2-D  station  was  used  to  calculate  the  skin- 
friction  coefficient  at  the  wall  using  Coles’  law-of-the-wall  fit,  therefore  the  profile  follows 

the  Coles’  law  of  the  wall  ( U*  =  — In  ( y*)  +5)  (Fig.  4).  In  tunnel  coordinates,  at  Station 

0.41 

5  the  U/tv^  profile  is  lower  above  y*=50.  At  the  separation  station  velocity  gradients  are 
higher  near  the  wall.  The  profile  is  also  higher  than  both  the  2D  and  Station  5  profiles. 
At  the  Vortex-core  station  in  the  region  y^=50  to  200  the  U/u,  velocity  gradient  is 
approximately  zero. 

The  V/u^  velocities  are  not  affected  by  the  coordinate  system  rotation  (Fig.  5).  At 
the  2-D  station  velocities  increase  towards  the  layer  edge.  At  Station  5  a  peak  at  y"=500 
forms.  At  separation  station  a  negative  peak  at  y*~60  is  accompanied  by  a  positive  peak 
at  y^*500.  At  the  Vortex-core  station  velocities  increase  till  y*=400  and  a  sharp  decrease 
towards  the  layer  edge  is  observed.  V/u^  changes  sign  and  the  magnitude  is  increased 
towards  the  Vortex-core  station. 

The  W/u,  at  2-D  station  is  very  close  to  zero  (Fig.  6 ).  The  velocity  magnitude  peak 
is  increased  and  shifted  further  away  from  the  wall  progressing  towards  the  Vortex-core 
station.  At  the  Vortex-core  station  the  W/u^  decreases  quickly  above  the  peak  location. 

I-B)  Normal  stresses 

At  Station  5  the  profile  values  are  1/4  lower  than  the  values  at  2-D  station 


throughout  the  layer  (Fig.  7 ).  However  at  the  Separation  station  the  values  are  increased 
to  follow  the  2-D  variation  very  closely  except  in  the  y''*=15  to  150  region.  At  the  Vortex- 
core  station  the  values  are  increased  by  1 .5  times  compared  to  2-D  and  about  2  times 

compared  to  the  Station  5  values.  Above  y"=400  a  quick  decrease  in  the  results 

in  values  less  than  the  2-D  values  above  y"=800. 

The  increases  to  very  large  values  at  the  Vortex-core  station  (Fig.  8  ). 

at  Station  5  compared  to  2-D  flow  near  wall  values  between  y"=20  to  150  region 

have  differences.  Between  y'^^OO  to  1000  region  values  at  Separation  station  are 

=1 .5  times  larger  than  the  2-D  values.  At  the  Vortex-core  station  is  ~7  times  larger 
than  2-D  value  peak. 

At  the  Vortex-core  station  the  values  are  also  different  from  the  2-D 

distribution  (Fig.  9  ).  The  peak  value  at  the  Vortex-core  station  is  8  to  9  times  larger  than 

the  2-D  values.  At  Station  5  near  the  wall  is  slightly  higher  than  the  2-D  values. 

At  the  Separation  station  values  are  2  to  3  times  to  that  of  2-D  values. 

The  profile  very  near  the  wall  and  above  y*«1000  vary  similarly  at  all 

stations.  The  TF/u/  show  similar  variations  above  y*=1000  at  every  station. 

I-C)  Shear  stresses 

The  Uvfu^^  shear  stress  profiles  in  the  3-D  flow  also  show  larger  variations  than 
the  2-D  flow  (Fig.  10).  At  Station  5  and  at  the  Separation  station,  near  wall  values  below 
y"«300  and  below  y*=1 00  are  lower  than  the  2-D  values,  respectively.  At  the  Separation 
station  above  y^’^lOO  the  values  are  higher  than  the  2-D.  Vortex  core  station  values  peak 
at  y'^=400  with  a  value  =3  times  higher  than  the  2-D  values. 

The  Uw  and  vw  shear  stresses  in  a  2-D  TBL  are  nominally  zero.  The  location  in 

the  3-D  flow  profiles  where  Uw/u^^  is  zero  shifts  towards  the  wall  progressing  towards  the 
Vortex-core  station  (Fig.  11).  The  positive  peak  value  at  the  Vortex-core  station  is  an 
order  of  magnitude  larger  than  the  Station  5  peak  value.  Above  y^-SOO  the  stress  values 

are  very  close  to  zero,  showing  that  the  TTvv/ u/  stress  in  this  region  is  not  effected  by  the 
near  wall  structure. 

Similarly,  viv  stresses  are  also  very  close  to  zero  above  y*~700  (Fig.  12).  The 

positive  and  the  negative  peak  values  are  increased  and  the  location  where  Vw~0  in  the 
profile  is  closer  to  the  wall  towards  the  Vortex-core  station. 

I-D)  Triple  products 

Triple  products  very  broadly  describe  the  flux  of  a  normal  stress  in  some  certain 

direction.  For  example  7Fi/  is  the  flux  of  TF  normal  stress  in  the  v  fluctuating  velocity 
direction.  Also  the  triple  products  give  information  about  the  simultaneous  occurrences 

of  different  velocity  fluctuations  in  the  average  sense.  For  example  7Fi/>0  shows  that  the 


u  fluctuating  velocities  in  the  average  are  associated  with  large  positive  v  fluctuating 

velocities.  Also  this  example  shows  that  at  the  locations  where  lFv>0  the  motion  is 
associated  with  positive  v  fluctuations,  therefore  moving  away  from  the  wall.  The  "sweep" 
and  "ejection"  motions  are  the  motions  when  (u>0,  v<0),  and  (u<0,  v>0)  fluctuation 

velocities  occur  simultaneously,  respectively.  Therefore  the  locations  where  'lFv<0  and 

uv^>0  denote  "sweep"  motion  dominated  regions.  Quadrants  defined  with  the  signs  of 
the  u  and  the  v  fluctuating  velocities  can  be  generalized  to  include  the  sign  of  w 
fluctuations  therefore  octants  can  be  defined. 

In  general  the  triple  products  at  the  Vortex-core  station  are  larger  compared  to  the 

other  stations  (Fig.  13).  The  variation  of  TFv  profiles  show  the  local  maxima  and  two  local 
minima  at  all  stations.  The  minimum  near  the  wall  is  located  at  y^=1 0  showing  that  the 

flux  of  IF  is  towards  the  wall.  Second  local  minimum  is  located  around  y^'^lOO  to  250 
region  shows  a  very  high  negative  value  at  the  Vortex-core  station.  At  the  Vortex-core 

station  large  gradients  of  the  quantity  are  observed.  Direction  of  the  IF  flux  changes  3 
times  at  this  station,  similar  to  that  at  Station  5  with  larger  magnitudes. 

The  IFw  is  nominally  zero  for  the  2-D  flow  (Fig.  14).  Large  negative  values  at 
y*=8  to  10  are  followed  by  large  positive  values  at  y^«235  to  300  at  all  stations.  Large 

differences  between  the  Separation  and  the  Vortex-core  station  are  observed  for  Ww/u^ . 
While  the  quantity  peaks  at  y^=200  at  the  Vortex-core  station  with  positive  values.  Station 

5  peak  values  are  negative  at  ^=300.  This  shows  that  the  lateral  flux  of  the  IF  normal 
stress  is  opposite  to  each  other  at  these  two  stations. 

The  is  also  nominally  zero  at  2-D  station  (Fig.  15).  However  the  values  at 
Separation  and  Vortex-core  stations  are  overwhelmingly  large  compared  to  Station  5. 

IFw/  shows  a  negative  peak  at  ^*250  at  Separation  station  where  it  shows  a  positive 

peak  at  y'’=45  and  a  negative  peak  at  y^=400  at  the  Vortex-core  station.  The  flux  of  IF 
with  the  lateral  fluctuations  change  sign  in  the  layer  at  y'=150  at  the  Vortex-core  station. 

The  IFu/u®  is  also  overwhelmingly  large  with  negative  values  at  the  Vortex-core 
station  compared  to  the  other  stations  (Fig.  16).  The  value  shows  a  negative  peak  at 
y‘"=400.  The  negative  peak  value  at  Separation  station  is  located  at  y"=225.  The  sign  of 
the  quantity  stays  the  same  through  out  the  Separation  and  Vortex-core  stations.Showing 

that  the  flux  of  IF  is  associated  with  the  negative  u  fluctuations.  Together  the  TFv  anduvF 
at  the  Vortex-core  station  show  that  near  the  wall  and  in  the  outer  region  ejection  motions 
are  highly  dominant  and  between  y^«150  to  300  negative  v  fluctuations  are  associated 
with  negative  u  fluctuations. 

The  IFF/  u®  values  show  different  characteristics  at  different  stations  (Fig.  1 7).  The 
location  where  the  quantity  changes  sign  from  positive  top  negative  near  the  wall  is 
closest  to  the  wall  proceeding  towards  the  Separation  station.  Values  above  y'^^lOO  at  2- 
D  and  Separation  stations  are  much  smaller  than  the  values  at  Separation  and  Vortex- 
core  stations.  At  the  Separation  station  and  at  the  Vortex-core  station  profiles  show  two 
negative  peaks.  At  the  Vortex-core  station  also  a  positive  peak  is  observed  at  y*«175. 

Magnitudes  of  the  FF/uf  values  are  very  small  compared  to  the  values  at 


separation  and  Vortex-core  stations  (Fig.  18).  Near  wail  values  at  3-D  stations  are 
negative  below  y*=25.  At  Separation  and  Vortex-core  stations  the  profiles  show  two 

peaks.  Negative  vw'^  shows  that  TP  normal  stress  flux  is  towards  the  wall.  Combined 

with  the  observations  for  Pw  figures  show  that  below  y*=25  there  is  a  region  where  v<0 
and  w>  0  fluctuations  are  dominant.  Between  y‘'«25  till  y'"=150,  w>0  is  accompanied  by 
v>0  and  between  y"^«1 50  till  y^==550  w<0  is  accompanied  with  v<0  fluctuations. 

The  uvw/u^  is  nominally  zero  at  the  2-D  station  (Fig.  19).  At  all  stations  near  wall 
values  below  y'^l  7  are  positive.  At  Separation  station  and  at  Station  5  the  values  change 
sign  to  positive  between  y'^^lOO  to  200.  However  the  variation  of  the  quantity  at  the 
Vortex-core  station  changes  sign  many  times  in  the  layer. 

The  TP/u^  profiles  show  positive  peaks  between  yW  to  15  at  different  stations 
(Fig.  20).  The  sign  of  the  variation  becomes  negative  at  different  y*  locations  between 
y''=15  to  30.  At  2-D  station  and  at  Station  5  the  quantity  becomes  positive  between  y'^^GO 
to  300.  However  values  at  the  Separation  station  and  the  Vortex-core  station  do  not 
change  sign,  stay  negative  and  show  two  local  minima. 

The  IP/u^  values  at  the  Vortex-core  station  are  overwhelmingly  larger  than  the  2- 
D  and  Station  5  values  (Fig,  21).  The  maximum  value  obtained  at  Vortex-core  station  is 
=14  times  larger  than  the  maximum  value  =1  obtained  at  Separation  station. 

The  TP/u®  values  reach  to  very  high  negative  values  (=-50)  around  y'"=30  at 
Vortex-core  station  (Fig.  22).  This  value  is  approximately  an  order  of  magnitude  larger 
than  the  ones  obtained  at  Separation  station  and  two  orders  of  magnitude  of  the  values 
at  Station  5.  Vortex-core  station  profiles  show  two  smaller  peaks  above  y‘"«1 00, 

1-E)  Estimates  of  uncertainty 

The  uncertainty  in  each  term  of  the  mean  velocity,  Reynolds  stress,  and  triple 
products  were  calculated  using  the  two  different  data  sets  acquired  at  the  Separation 
station  and  using  Chauvent’s  criterion  to  calculate  the  standard  deviation.  The  equation 
used  is 

^  =  1.15 

G 

where  d^^g^  'S  the  average  of  the  half  of  the  differences  between  two  calculated  values  for 
that  term.  20  to  1  odds  uncertainties  calculated  as  120  are  tabulated  in  table  1. 

I-F)  Derived  quantities 

l-F-1)  Flow  angie  (FA) 

The  flow  angle  is  calculated  using, 

FA  =arctan(— ) 

U 

The  2-D  flow  angles  are  nominally  zero  (Fig.  23).  Flow  angle  magnitudes  increase 
proceeding  towards  the  Vortex-core  station.  The  minimum  location  in  the  profile  shift  up 
into  the  layers.  Around  y"=800  flow  angles  at  every  station  are  close  to  one  another.  Flow 


angle  at  Vortex-core  station  increases  more  steeply  at  y"  above  the  minimum  location. 


l-F-2)  Flow  Gradient  angle  (FGA) 

The  FGA  is  calculated  using, 

FGA  =arctan(.|5^^) 

dU !  dy 

The  FGA  variation  near  the  wall  below  y^«10  is  similar  at  Station  5  and  at  Separation 
station  (Fig.  24).  The  values  at  the  Vortex-core  station  near  the  wall  are  lower  than  the 
other  two  stations.  The  FGA  changes  sign  at  higher  y*  locations  in  the  profiles.  FGA  at 
the  Vortex-core  station  passes  through  zero  value  vary  abruptly.  The  magnitude  of  the 
values  are  larger  proceeding  towards  the  Vortex-core  station. 

l-F-3)  Shear  stress  angle  (SSA) 

The  shear  stress  angle  is  calculated  using, 

SSA  =arctan  ( _ ^ ) 

-uv 

The  magnitude  of  the  SSA  in  the  profiles  increases  and  the  locations  in  the  profiles  where 
SSA  changes  sign  are  closer  to  the  wall  proceeding  towards  the  Vortex-core  station  (Fig. 
25).  The  FGA  and  SSA  profile  values  do  not  follow  each  other.  While  at  Station  5,  FGA 
and  SSA  change  sign  at  a  location  closer  to  each  other  at  the  Separation  station  and  at 
Vortex-core  station  SSA  changes  sign  closer  to  the  wall.  Near  wall  values  of  SSA  are 
much  different  than  the  FGA  values. 

l-F-4)  A^,  Townsend’s  structural  parameter 

This  structural  parameter  is  defined  as; 

IF  +'v^  +T^ 

The  variations  of  the  parameter  below  y'^IO  are  close  to  each  other  at  the  2-D,  Station 
5  and  Separation  stations  (Fig.  26).  At  the  Vortex-core  station  this  variation  is  not 
followed.  Between  y«100  to  1000  the  value  is  close  to  a  constant  at  the  2-D  flow  station. 
At  Station  5  the  values  reach  to  a  maximum  of  0.17  at  y^=1 000.  At  the  Separation  station 
the  values  are  reduced  throughout  the  layer  compared  to  Station  5.  At  the  Vortex-core 
station  very  large  double  peaks  are  observed.  This  shows  that  u,  v,  and  v,  w  fluctuating 
velocities  are  better  correlated  for  this  station. 

l-F-5)  Turbulent  Kinetic  Energy  (TKE) 

The  TKE/u/  is  calculated  using 

T/CE^IF+'iF+TiF 

The  TKE  of  the  flow  at  Station  5  is  reduced  compared  to  2-D  station  (Fig.  27).  The  values 


rapidly  increase  towards  the  Vortex-core  station.  At  the  Vortex-core  station  values  are 
=2.5  times  that  of  the  2-D  station.  Near  the  wall  and  above  ^=1000,  variations  are  similar 
to  each  other  at  every  station. 

l-F-6)  Anisotropy  constant  (Njc) 

The  anisotropy  constant  Njc  is  calculated  using, 

^-vw/{dWtdy) 

-uvlidUldy) 

The  Ntc  is  the  ratio  of  the  eddy  viscosities  in  the  z  and  the  x  directions  (Fig.  28).  Many 
algebraic-eddy  viscosity  turbulence  models  assume  that  the  ratio  is  =1 .  The  figure  shows 
that  the  ratio  is  increasingly  away  from  1  towards  the  Vortex-core  station  near  the  wall. 
Towards  the  layer  outer  edge  the  values  become  close  to  zero. 


CHAPTER  II  -  TRANSPORT  BALANCE  OF  THE  STRESSES 


The  Reynolds-stress  transport  equation  with  summation-subscript  notation  can  be 
written  as: 

where 


"  3t 


convection 


, _ dU:  _ 94, 


n,dU,  dU:^ 

<D..=i-( _ 1+ _ i) 

“  p^9x^.  dx' 

d(Muiu,) 

dx, 

“  P  Sx,  -- 

d'^U.U, 


TD,r- 


dX; 


VD„  =  ^v—L^ 
‘I  .2 


dx; 


production 
pressure-strain 
turbulent  diffusion 
)  pressure-diffusion 
viscous-diffusion 


^  du,  du, 

e..  =  +2v _ - — i 

“  dx,  dx. 


viscous-dissipation 


The  subscripts  denote  the  axes.  Repeated  dummy  indices  in  a  multiplication  denote 
summation.  Here  UJJj  denote  the  Reynolds  stress  tensor,  Uj  stands  for  the  three 
components  of  the  mean  velocity  vector,  p  is  the  fluctuating  pressure,  v  is  the  kinematic 
viscosity,  p  is  the  density,  and  t  is  the  time.  Note  the  minus  sign  in  front  of  the  e,j. 

All  the  terms  are  point  dependent  except  the  and  PD,j  which  are  functions  of 
space  because  of  the  pressure  fluctuation.  The  terms  with  the  fluctuating  pressure  contain 
the  history  effects  due  to  the  mathematical  nature  of  the  pressure  fluctuations.  Pressure 
fluctuations  at  a  point  can  be  expressed  as  the  integral  of  correlations  of  velocity 
fluctuations  in  the  whole  domain  using  the  Green’s  theorem  on  solving  the  Poisson 
equation  (Rotta,  1962).  Measurement  techniques  measurement  of  pressure  fluctuations 
within  a  flow  is  not  possible  at  high  Reynolds  numbers.  Measurement  of  viscous 
dissipation  with  reasonable  accuracy  very  near  a  boundary  proves  to  be  very  difficult. 
Therefore  in  the  individual  stress  transport  equations  PD,j,  e,j  terms  can  not  be 
directly  measured  or  calculated  from  the  measured  data.  While  approximate  approaches 
are  to  be  used  for  two  of  the  terms,  the  third  one  can  be  extracted  from  the  balance  of 

the  transport  equation.  The  turbulent  kinetic  energy  (TKE=.lu,  u,)  transport  equation  is 

equal  to  the  sum  of  the  transport  equations  for  the  normal  stresses  divided  by  2 
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(Convection  of  TKE= 


C33)).  Due  to  the  continuity  equation  for  incompressible 


flows  ( — i=0),  the  pressure-strain  term  in  the  TKE  transport  equals  zero.  Therefore  TKE 

transport  equation  involves  two  unknown  terms;  pressure-diffusion  and  viscous 
dissipation.  This  allows  one  to  use  an  approximate  model  for  one  of  the  terms  and  extract 
the  other  term  from  the  balance  of  the  TKE  transport.  Generally  the  pressure-diffusion 
term  is  accepted  to  be  very  small  throughout  the  layer  and  an  error  is  therefore 
introduced  to  the  dissipation  term  (Schwarz  and  Bradshaw,  1994). 

The  pressure-diffusion  term  for  a  homogenous  flow  was  modelled  by  Lumley 
(1978)  using  a  Gaussian  model.  In  this  model  he  made  use  of  the  moment  generating 
function  to  generate  an  equivalent  scalar,  velocity  joint  density.  Using  this  equation  and 
using  Fourier  transforms  for  homogeneous  stochastic  fields  (Lumley,  1970),  he  obtained 
the  equation 

P  5 

However,  the  large-eddy  simulations  of  Shao  et  al.  (1990),  of  a  shearless  turbulence 
mixing  layer  flow  shows  that  the  constant  (1/5)  is  an  overestimated  figure,  and  neglecting 
pressure  diffusion  gives  a  better  prediction.  The  constant  multiplier  was  also  discussed 
in  Shih  et  al.  (1987),  but  a  new  number  was  not  proposed.  The  pressure-diffusion  term 
was  assumed  to  be  small  by  Schwarz  and  Bradshaw  (1993,  1994)  mainly  because  their 
work  compared  the  experimentally  extracted  pressure-strain  data  to  the  theoretical  ones 
above  y'^^ISO.  Indeed  the  pressure-diffusion  terms  prove  to  be  small  away  from  the  wall 
above  y*=100.  However  the  pressure-diffusion  term  is  not  necessarily  small  below  this  y* 
location  using  Lumley’s  approximation.  In  fact  the  study  of  Launder  and  Tselepidakis 
(1990)  of  second  moment  closure  within  the  sublayer  show  that  the  pressure  diffusion  has 
substantial  importance  in  shaping  the  sublayer  turbulence  structure.  They  used  0.15  for 
the  constant  coefficient  but  the  constant  was  not  derived  through  a  rigorous  analysis. 
Latest  DNS  studies  (Sumitani,  and  Kasagi,  1995)  show  that  the  pressure-diffusion  term 
in  the  normal  and  shear  stresses  are  large  and  as  large  as  the  pressure-strain  term  below 
y'^IO  near  the  wall  and  the  pressure-  diffusion  term  for  the  TKE  balance  is  practically 
zero.  This  is  also  approximately  true  for  the  pressure-diffusion  terms  estimated  for  the 
TKE  equation  using  Lumley’s  approximation.  This  ensures  that  the  dissipation  term  in  the 


TKE  transport  equation  £=v — i  — -  calculated  whether  using  Lumley’s  approach  or  not 

3x,  dx, 

will  not  be  very  much  in  error.  Since  the  pressure-strain  models  which  will  be  discussed 
later  mainly  use  the  e  and  the  anisotropy  of  the  shear  stress  tensor,  they  are  also  not 
affected  by  this  approximation  to  the  pressure-diffusion  term.  However,  the  individual 
pressure-strain  terms  extracted  from  the  data  will  have  large  errors  especially  near  the 
wall  if  the  pressure-diffusion  terms  are  neglected.  Therefore,  the  error  introduced  by  using 
this  approximation  is  believed  to  be  smaller  than  neglecting  it.  If  the  term  is  neglected, 
the  O;.  calculated  will  include  the  effect  of  the  pressure-diffusion  term  in  it. 

The  TKE  dissipation  was  extracted  from  the  TKE  transport  equation  budget 
balance.  In  order  to  calculate  the  pressure  diffusion  both  in  TKE  transport  equation  as 
well  as  in  the  stress  transport  equations,  Lumley’s  approximation  with  its  original 
coefficient  was  used.  To  judge  the  effect  of  neglecting  or  keeping  the  pressure-diffusion 


term  in  the  balance  of  the  TKE  budget,  the  very  near  wall  variation  of  the  dissipation  term 
could  not  be  extracted  from  the  measured  normal  and  shear  stress  data.  The  near  wall 
limiting  value  of  the  dissipation  and  variation  very  near  the  wall  can  be  written  using  the 
continuity  relations  and  the  variations  of  the  fluctuating  velocities  very  near  the  wall.  The 
analysis  showed  that  in  the  calculations  of  the  near  wall  variation  of  the  dissipation,  at 
least  second  order  approximations  perpendicular  to  the  wall  need  to  be  made  in  order  for 
the  variation  to  be  physically  possible.  This  second  order  approximation  could  not  be 
calculated  from  the  measured  normal  and  shear  stress  data  since  the  coefficients 
extracted  from  near  wall  variations  of  normal  and  shear  stresses  used  in  the  dissipation 
variation  become  coupled.  The  analysis  of  E/j  near  the  wall  shows  that  (Reynolds,  1 983, 

Lai,  and  So,  1990,  Mansour,  et  al.,  1988)  the  dissipation  of  the  IF  and  wF  vary  linearly 

with  y  near  the  wall  with  a  constant  value  at  the  wall  and  that  the  dissipation  for  thev^ 
is  zero  at  the  wall  and  varies  as  y^  near  the  wall.  This  necessitates  use  of  an  anisotropic 
dissipation  tensor  for  asymptotically  correct  modelling.  However  the  correction  factor  for 
the  present  study  was  calculated  as  practically  zero  using  the  anisotropic  dissipation 
tensor  of  Lai  and  So,  (1990).  In  this  study  the  effect  of  anisotropic  dissipation 
approximation  on  the  stress  transport  budget  and  on  the  predictions  of  the  pressure-strain 
models  is  also  tested  and  reported  using  the  anisotropic  dissipation  rate  algebraic  model 
of  Hallback,  Groth,  and  Johansson  (1990). 

Viscous  dissipation  terms  in  shear  stress  transport  equations  are  found  to  be  small 
by  direct  numerical  studies  of  3-D  TBLs  at  low  Reynolds  numbers  (Spalart,  1989). 
However,  the  viscous  dissipation  terms  in  the  normal  stress  transport  equations  are  not 
small. 

Once  an  isotropic  dissipation  is  assumed,  the  original  pressure-diffusion  terms 
included  in  the  dissipation  terms  for  the  individual  stresses  are  replaced  with  equal 
distribution  for  each  of  the  stresses.  Use  of  an  anisotropic  dissipation  tensor  is  still  a 
research  topic.  At  high  Reynolds  numbers  the  dissipation  is  usually  assumed  isotropic. 
However  within  the  limits  of  the  proposed  algebraic  dissipation  rate  tensors  (Hallback, 
1994),  an  anisotropic  dissipation  approximation  is  much  different  than  the  isotropic 
dissipation.  Under  the  restrictions  of  present  discussions  it  is  believed  that  the  pressure- 
diffusion  term  should  be  kept  at  least  as  proposed  by  Lumley. 

The  dissipation  in  this  study  was  estimated  both  using  the  isotropic  and  anisotropic 
dissipation  approximations.  The  models  which  account  for  the  anisotropy  of  the 
dissipation,  show  effects  only  very  near  the  wall  below  y'^=5  due  to  very  small  correction 
factor  f^i  as  proposed  by  Lai  and  So  as  discussed  in  Chapter  V,  or  to  very  large 
differences  in  the  dissipation  rates  throughout  the  layers  as  proposed  by  Hallback  et  al.. 
7he  effect  of  anisotropic  dissipation  rate  and  near  wall  corrections  are  discussed  later  in 
Chapter  IV. 

From  an  experimental  point  of  view  derivatives  of  any  quantity  can  not  be 
measured  directly,  therefore  approximations  are  made  using  either  a  Taylor  series 
expansion  or  fitting  a  function  to  the  data  and  calculating  the  derivatives  of  this  function 
at  that  point  of  interest.  The  gradients  of  the  mean  velocities  and  stresses  in  the  y 
direction  were  calculated  using  a  parabola  fit  at  5  consecutive  points  to  the  original  data 
and  taking  the  derivative  of  this  curve  fit  at  the  third  point  location.  As  mentioned  earlier, 
several  profiles  at  closely  spaced  x  and  z  separations  around  the  four  main  measurement 
locations  were  made.  Derivatives  in  the  y  direction  were  calculated  at  each  measurement 
location  separately.  The  original  data  values  were  next  replaced  with  the  parabola  fit 
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values  after  the  derivatives  of  the  1 9  quantities  (mean  velocities,  stresses,  and  the  triple 
products)  were  calculated  in  the  y  direction.  The  procedure  was  used  in  order  to  smooth 
the  data  and  therefore  to  reduce  the  uncertainties.  Next,  y  locations  of  one  of  the  profiles 
were  chosen  as  the  base.  Other  profile  data  values  were  replaced  with  linear  interpolation 
results  at  the  same  y  locations  of  the  base  profile.  The  x  and  2  derivatives  were 
calculated  either  with  a  second  or  a  first  order  accurate  in  space  finite  difference 
approximation  using  the  interpolated  data. 

Calculations  of  the  "Convection",  "Production",  "Turbulent  diffusion",  and  the 
"Viscous  diffusion"  terms  are  straightforward.  The  Viscous-diffusion  term  was  calculated 
best  by  a  5  point  parabola  fit  to  the  first  order  y  derivative  of  the  variables  and  calculating 
the  derivative  of  this  parabola  fit,  rather  than  a  5  point  parabola  fit  to  the  data  values  and 
taking  the  second  derivative  of  this  function.  Another  scheme  tested  was  a  5  point 
forward-finite-differencing  scheme  which  resulted  in  high  fluctuations  of  the  second  order 
derivatives.  In  the  calculation  of  the  viscous  diffusion  term  only  the  y  derivatives  were 
used  to  reduce  the  effect  of  the  uncertainties  involved  in  calculating  the  second  order 
derivatives  in  x  and  z  directions.  The  first  order  derivatives  used  are  the  derivatives 
interpolated  to  the  same  y  locations  of  the  base  stations. 

The  £  term  in  the  TKE  equation  was  extracted  using  Lumley’s  approach  for  the 
pressure-diffusion  term.  The  "Dissipation"  terms  in  the  normal  stress  transport  equations 
were  approximated  both  with  isotropic  dissipation  and  the  anisotropic  dissipation  model 
of  Hallback,  Groth  and  Johansson  (1991).  In  the  isotropic  dissipation  assumption  for 
individual  normal  stresses  2/3  of  the  TKE  Dissipation  was  used.  The  "Pressure-Strain" 
terms  in  the  both  normal  and  shear  stress  equations  were  extracted  from  the  balance  of 
the  equations.  The  dVIdy  term  was  extracted  from  the  continuity  equation  for  the  mean 
flow. 

Figures  28  through  52  show  the  terms  in  the  transport  of  individual  stresses 

nondimensionalized  with  u^N  vs.  in  the  tunnel  coordinates.  For  each  flow  type,  the  six 
plots  presented  belong  to  the  terms  in  the  transport  equation  for  the  Reynolds’  stresses. 
In  this  section  the  dissipation  is  assumed  isotropic,  and  therefore,  equally  distributed 
among  the  normal  stresses.  Dissipation  terms  in  the  shear  stress  equations  are  assumed 
zero.  Individual  terms  with  anisotropic  dissipation  approximation  are  presented  and 
discussed  in  Chapter  IV.  Pressure-diffusion  terms  are  estimated  using  Lumley’s  approach. 

Il-A)  Estimates  of  uncertainty 

The  uncertainty  in  each  term  of  the  stress  balance  equations  were  calculated  using 
the  two  values  calculated  for  that  individual  term  using  the  two  different  data  sets 
acquired  at  the  Separation  station  and  using  Chauvent’s  criterion  to  calculate  the  standard 

deviation.  20  to  1  odds  uncertainties  calculated  as  ±2  a  are  tabulated  in  table  2. 

ii-B)  Stress  transport  equation  budget 
ii-B-1)  Baiance  at  the  2-D  station 

Transport  equation  budget  for  IF  stress  shows  that  the  production  term  is  the 
largest  term  among  all  terms  especially  near  the  wall  region  (Fig.  29).  The  IF  production 
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peaks  near  the  wall  at  y^=15  to  20  with  a  value  of  0.34.  The  turbulent  diffusion  changes 
sign  at  y^=40  and  magnitude  at  the  peak  location  (y*«20)  is  comparable  to  dissipation  and 
convection  peak  magnitude  terms.  Above  y'^^tOO  pressure-strain,  production  and  the 
dissipation  terms  are  larger  than  other  terms.  The  pressure-diffusion  term  is  practically 
zero  throughout  the  layer  for  this  stress.  Except  the  convection  term  ail  profiles  show 
peaks  between  y'^^t  5  to  20  range. 

The  transport  balance  shows  that  the  major  contributors  are  the  pressure-strain 
and  the  dissipation  terms  with  the  same  sign  (Fig.  30).  The  F  production  is  practically 
zero  for  the  2-D  station  since  the  major  contributors  -UvdVIdy  and  -'^dVIdy  have 

opposite  signs  to  each  other  throughout  the  layer  and  v^/  |t/p  j  =1  for  2-D  boundary 
layers.  Turbulent  diffusion  changes  sign  around  y^«=10,  20,  and  100.  Convection  value 
peaks  around  y^=15  and  both  are  important  only  near  the  wall  below  y^=100.  Above 
y*=100  dissipation  and  pressure-strain  terms  follow  each  other  closely.  In  the  near-wall 
region  below  y*«50  pressure-diffusion  term  is  a  large  contributor  with  a  peak  around 
y*=20  with  a  value  =2/5th  of  the  dissipation  term.  Viscous  diffusion  is  important  only  below 
y^=20. 

Production  of  is  practically  zero  in  a  2-D  flow  since  1jw,vw,  and  dWldy  are 
zero  (Fig.  31).  The  pressure-strain  term  is  the  highest  term  in  the  transport  of  TF.  This 

also  shows  that  the  TF  normal  stress  is  due  to  redistribution  of  the  TKE  by  the  pressure- 
strain  term  between  the  v  and  w  fluctuations.  Convection,  viscous  diffusion,  and  turbulent 
diffusion  terms  are  important  near  the  wall  below  y''=40.  Dissipation  values  and  pressure- 
strain  values  follow  each  other  closely  above  y^=50. 

For  the  normal  stresses  in  general  the  dissipation  and  the  pressure-strain  terms 

are  at  the  same  order  of  magnitude.  Pressure-strain  term  redistributes  the  7F  production 
to  the  other  fluctuating  components. 

In  this  section  of  the  discussion  the  dissipation  in  Uv  transport  equation  is 
approximated  by  zero,  based  on  the  DNS  simulations.  The  production,  pressure-strain 
and  the  pressure-diffusion  terms  are  the  dominating  terms  throughout  the  layer.  An  order 
of  magnitude  smaller  convection,  turbulent-diffusion,  viscous-diffusion  terms  are  away 
from  zero  below  ^=50  to  60  (Fig.  32).  Pressure-diffusion  term  peaks  at  the  same  location 
of  the  production  term  at  y'^^IS  with  a  value  half  of  that  of  the  production  term.  Other 
stresses  in  a  2-D  flow  are  zero  (Figures  33,  and  34). 

Turbulent  kinetic  energy  budget  shows  that  the  dissipation  and  the  production 
terms  follow  each  other  closely  above  y^«35,  and  all  other  terms  are  close  to  zero  above 
this  location  (Fig.  35).  Production  profile  peaks  at  y*=15  to  20  and  values  are  larger  than 
the  dissipation  term  values  at  this  location.  Pressure-diffusion  and  viscous-diffusion  have 
similar  profiles  with  opposite  signs  to  each  other.  Compared  to  the  production  and 
dissipation  terms  the  pressure-diffusion  is  much  smaller  through  out  the  layer,  which 
shows  that  the  dissipation  term  calculated  from  the  TKE  budget  balance  using  Lumley’s 
approach  for  the  pressure-diffusion  term  does  not  have  a  substantial  effect  on  this  value. 
The  convection  term  at  y''»10  is  close  to  the  dissipation  term  values.  Turbulent  diffusion 
values  are  twice  that  of  the  viscous  diffusion  at  y^==20. 

ll-B-2)  Balance  at  Station  5 
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At  Station  5  TP  transport  budget  shows  that  production  term  is  higher  than  the 
other  terms  (Fig.  36).  All  the  term  values  are  close  to  zero  above  y'"=125.  Below  this 
location  turbulent-diffusion,  dissipation,  convection,  and  the  pressure-strain  are  all  large 
and  important.  Viscous-diffusion  term  is  important  below  y''«30  are  on  the  same  order  of 
magnitude.  Compared  to  the  2-D  flow  all  the  terms  are  reduced  in  magnitude.  The 
decrease  in  magnitude  for  pressure-strain  and  production  terms  are  more  pronounced 
than  the  other  terms. 

In  the  transport  of  TP,  pressure-strain,  dissipation  and  pressure-diffusion  terms  are 
the  largest  contributors  (Fig.  37).  Compared  to  2-D  flow  the  convection  and  turbulent 
diffusion  terms  are  much  smaller.  Pressure-diffusion  term  shows  a  peak  at  y'=20  and  is 
negligibly  small  above  y"^=100.  Above  y''=100  the  dissipation  and  pressure-strain  term 
values  are  close  in  magnitude. 

On  contrary  to  2-D  flow,  at  Station  5  production  of  TP  is  the  largest  term  in  the 
transport  budget  of  this  stress  (Fig.  38).  Dissipation  values  are  at  the  same  order  of 
magnitude  of  the  production  but  are  smaller  throughout  the  layer.  Convection,  pressure- 
strain,  turbulent-diffusion  and  viscous-diffusion  terms  are  =4  times  smaller  than  the 
production  throughout  the  layer.  Pressure-strain  and  turbulent-diffusion  term  profiles  vary 
as  mirror  images  to  one  another. 

The  pressure-diffusion,  pressure-strain,  and  production  terms  are  the  major  terms 

in  the  budget  balance  of  the  Uv  stress,  and  the  turbulent  diffusion,  convection  and  the 
viscous-diffusion  are  negligibly  small  (Fig.  39).  Pressure-strain  term  variation  follows  the 
pressure-diffusion  term  variation  with  an  opposite  sign. 

The  Uw  production  is  larger  than  the  Uv  production  (Fig.  40).  Pressure-strain  and 
production  terms  are  the  major  terms  and  they  have  opposite  signs  throughout  the  layer. 
The  turbulent-diffusion,  viscous-diffusion  and  convection  terms  are  larger  in  magnitude 

then  the  equivalent  Uv  terms. 

The  convection  and  turbulent  diffusion  terms  are  also  small  for  the  vw  transport 
equation  (Fig.  41).  The  dominant  terms  are  the  pressure-strain,  production,  and  the 
pressure-diffusion  terms.  The  distribution  of  pressure-diffusion  term  largely  shapes  the 
variation  of  the  pressure-strain  term. 

The  dissipation  and  the  production  terms  are  the  major  terms  in  the  balance  of  the 
TKE  (Fig.  42).  Production  and  dissipation  profiles  peak  at  y'^IS.  Effects  of  convection, 
turbulent,  viscous,  and  pressure-diffusion  terms  are  small  but  not  negligible  below  y"«40. 
The  diffusion  terms  change  sign  very  near  the  wall  in  the  y'^^S  to  12  range  and  also  at 
y+=40.  The  pressure-diffusion  and  viscous  diffusion  profile  peak  values  are  close  to  each 
other  with  opposite  signs. 

il-B-3)  Balance  at  Separation  station 

The  production  term  for  the  U^  balance  at  the  Separation  station  is  higher  than  the 
2-D  and  the  Station  5  values  near  the  wall  (Fig.  43).  The  dissipation  and  turbulent 
diffusion  are  increased  more  than  twice.  The  convection  term  does  not  change  sign  within 
the  layer  as  observed  in  2-D  station  and  is  positive.  However  all  the  terms  reduce  close 
to  zero  values  much  earlier  than  the  2-D  and  the  Station  5  values  at  y'^^SO.  Pressure 
diffusion  is  negligibly  small  and  viscous  diffusion  is  increased  especially  very  near  the 


OMADTPD  II 


13 


wall.  Large  pressure-strain  contribution  is  confined  into  a  narrow  region  near  the  wall  to 
y*=20  range  with  values  slightly  higher  than  that  of  Station  5. 

In  the  transport  equation  of  ^  the  dissipation  and  the  pressure-strain  terms  are 
increased  significantly  by  twice  compared  to  both  2-D  station  and  Station  5  values  (Fig. 
44).  The  production  is  also  increased  slightly  in  the  outer  region  around  ^=200.  The 
convection,  and  turbulent-diffusion  values  are  close  to  zero  at  this  station  and  are  slightly 
larger  than  zero  at  y*=300.  The  pressure-diffusion  term  is  =1/3  of  the  dissipation  term  at 
y*=20  with  larger  values  than  at  both  2-D  and  Station  5  profiles. 

Except  the  pressure-strain  the  terms  all  the  terms  in  the  W  budget  are  increased 
by  2.5  times  that  of  Station  5  (Fig.  45).  Pressure-strain  term  is  increased  by  5  times 
throughout  the  layer.  The  dissipation  and  the  production  terms  peak  at  y'=15  to  20.  The 
convection  term  is  negative  in  the  y*~7  to  300  range.  Pressure-strain  term  is  also 
negative  throughout  the  layer.  While  at  the  2-D  station  the  pressure-strain  and  dissipation 

are  the  main  terms  in  the  development  of  at  3-D  stations  the  production  plays  an 
important  role  diminishing  the  effect  of  the  pressure-strain. 

Compared  to  the  2-D  flow  the  pressure-strain  term  of  Uv  budget  is  doubled  and 
the  pressure-diffusion  term  is  tripled  (Fig.  46).  Production,  pressure-strain  and  pressure- 
diffusion  terms  all  show  extremum  at  y^=15.  The  location  where  the  pressure-diffusion 
term  changes  sign  very  near  the  wall  is  shifted  towards  the  wall.  With  increased  3-D 
effects  the  large  variations  near  the  wall  are  confined  to  a  smaller  near  wall  region. 
Convection,  turbulent-diffusion,  and  viscous-diffusion  terms  are  relatively  much  small 
throughout  the  layers. 

The  production  and  the  pressure-strain  terms  in  the  Uw  transport  equation  are 
increased  more  than  twice  compared  to  the  Station  5  values  at  y*=15  where  they  peak 
(Fig.  47).  The  turbulent  diffusion  is  also  increased,  especially  below  y^=10.  This  term  is 
large  with  values  as  high  as  1/5  of  the  production  peak  value  below  y*=40.  Viscous 
diffusion  increases  when  the  wall  is  approached  below  ^=30. 

Similar  observation  are  made  in  the  transport  equation  for  the  vw  as  in  the 

transport  equation  of  uv  (Fig.  48).  The  magnitudes  of  the  pressure-strain,  pressure- 
diffusion  and  the  production  terms  are  increased  by  more  than  twice  below  y^=100  and 
about  5  times  above  this  y*  location  compared  to  the  Station  5  values.  Pressure-diffusion 
term  below  y*=30  is  a  major  contributor  shaping  the  variation  of  the  pressure-strain 
profile.  The  production  and  the  pressure-strain  terms  change  sign  at  y"=100. 

TKE  budget  balance  shows  that  except  the  convection  term  all  the  term  values  are 
elevated  more  than  twice  compared  to  those  at  Station  5  (Fig.  49).  Effects  of  the  diffusion 
terms  are  important  below  y^=40.  Turbulent  and  pressure-diffusion  terms  change  sign  at 
y^=10  and  away  from  the  wall  at  ^=40.  The  profile  shapes  at  this  Separation  station  are 
not  highly  different  than  profiles  at  Station  5. 

ll-B-4)  Balance  at  Vortex-core  station 

The  peak  value  of  the  production  term  of  TP  shows  a  decrease  at  the  Vortex-core 
station  compared  both  to  the  2-D  and  the  other  3-D  cases  (Fig.  50).  The  peak  value  is 
comparable  to  the  Separation  station  value.  The  term  changes  sign  and  becomes 
negative  at  y*=60.  The  turbulent  diffusion  term  magnitude  is  as  high  as  the  value  at 
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Separation  station.  The  term  changes  sign  to  positive  at  y"=10  similar  to  the  other  station 
profiles.  The  pressure-strain  term  is  positive  near  the  wall  till  ^*125  on  contrary  to  the 
variations  at  other  stations.  The  dissipation  and  the  convection  terms  are  larger  in 
magnitude  compared  to  that  of  the  Separation  station.  Convection  profile  forms  a  negative 
peak  at  y''=30.  The  magnitudes  of  all  the  terms  are  not  small  in  the  outer  region  contrary 
to  the  observations  made  at  other  stations.  Viscous-diffusion  term  is  on  the  same  order 
of  magnitude  of  the  Separation  station.  Pressure  diffusion  is  different  than  zero  at  this 
station. 

The  profiles  of  the  terms  in  the  transport  of  are  highly  different  than  the  other 
station  profiles  (Fig.  51).  The  dissipation  term  shows  peaks  around  y*=30  and  around 
y*«300.  The  pressure-strain  term  is  approximately  flat  till  y'*’=150.  The  convection  term 
oscillates  around  zero.  The  turbulent  diffusion  term  peak  around  y^=150  and  magnitude 
of  the  peak  value  is  larger  then  that  of  the  production.  The  production  term  is  negative 
near  the  wall  and  changes  sign  around  y‘'=300  and  shows  a  positive  peak  around  y''=400. 
Pressure-diffusion  term  peak  value  at  y'"=25  is  increased  =2  times  that  of  Separation 
station. 

The  transport  terms  except  the  dissipation  term  show  a  further  increase  in 
magnitude  in  the  whole  layer  compared  to  the  Separation  station  (Fig.  52).  The 
dissipation  term  values  are  close  to  that  of  the  Separation  station.  The  elevated  levels  of 
the  terms  persist  further  up  in  to  the  layer.  Both  the  production  and  pressure-strain  values 
are  increased  more  than  twice  compared  to  the  Separation  station.  The  turbulent  diffusion 
term  attains  values  larger  than  the  pressure-strain  term  in  magnitude  at  y^=35  to  100 
region.  The  convection  term  shows  a  peak  minimum  at  y"=45  and  the  values  in  the  layer 
are  approximately  an  order  of  magnitude  larger  than  the  Station  5  and  2-D  station  values. 
Pressure-diffusion  is  positive  and  large  between  y^=20  to  300.  Viscous  diffusion  increases 
below  y"^®  10. 

The  terms  in  the  transport  equation  of  Uv  are  similar  in  shape  of  other  stations  with 
elevated  values  till  y''=50  (Fig.  53).  The  large  vortex  in  the  outer  region  effects  the 
variations  at  this  region.  Pressure-strain,  pressure-diffusion  and  production  terms  change 
sign  3  times  with  in  the  layer.  Production  and  pressure-strain  terms  show  local  peaks  with 
opposite  signs  at  y*~20  and  at  y''«400.  Pressure-diffusion  peak  magnitude  value  at  y*~2Q 
is  more  than  twice  of  the  production  peak  value  at  the  same  location. 

Near  wall  variations  of  the  Uw  transport  terms  below  y*=70  is  similar  to  other  3-D 
stations  (Fig.  54).  The  turbulent  diffusion  term  is  suppressed.  Near  the  vortex  core  at 
y'^=300  the  production  and  the  pressure-strain  values  show  large  local  peaks  with 
opposite  signs. 

One  of  the  largest  effects  of  the  vortex  located  in  the  outer  region  is  observed  in  viv 
transport  (Fig.  55).  The  production  term  distribution  term  is  similar  to  Station  5  values  with 
larger  negative  values  in  the  outer  region.  Between  to  150  pressure-diffusion  and 
pressure-strain  terms  vary  similarly  with  opposite  signs  to  each  other.  A  similar  variation 
between  pressure-strain  and  production  is  observed  between  y'"«150  to  800.  Pressure- 
strain  and  pressure-diffusion  peak  locations  for  Separation  station  at  y''=15  to  20  are 
shifted  away  from  the  wall  to  y*«60.  Near  the  wall  new  peaks  at  y'^^IO  are  developed. 
Pressure-strain  profile  shows  a  wave-like  variation  throughout  the  layer. 

The  profiles  of  the  terms  in  the  TKE  balance  at  this  location  differ  highly  than  the 
other  three  cases  studied  (Fig.  56).  Even  though  the  production  and  dissipation  values 


are  close  to  the  values  at  Separation  station  the  convection,  turbulent  diffusion,  pressure- 
diffusion  terms  are  highly  increased.  All  the  terms  except  the  viscous  diffusion  are 
important  in  the  whole  layers,  and  the  effect  of  the  terms  are  not  confined  to  y""  <  100 
region.  Turbulent  diffusion  terms  above  y^=1 00  are  comparable  to  the  production  and 
dissipation  terms.  Viscous  diffusion  changes  sign  at  y'^^l  0  and  increases  towards  the  wall. 
It  is  important  below  y'^OO.  High  negative  convection  values  are  obtained  away  from  the 
wall  with  a  peak  at  y*=40. 

Il-C)  Differences  in  the  stress  transport  due  to  fiow  conditions 
at  different  stations 

IF 

In  general  three  dimensionality  of  the  flow  reduces  the  production  of  the  IF  stress 
for  highly  skewed  3DTBLs.  At  Separation  station  however  the  production  is  higher  than 
the  2-D  flows.  Turbulent  diffusion  decreases  at  Station  5  compared  to  2-D  flow  station  but 
at  Separation  and  Vortex-core  stations  it  is  larger.  Turbulent  diffusion  progressively 
increases  towards  the  Vortex-core  from  Station  5.  Dissipation  terms  are  unchanged  at 
Station  5  compared  to  2-D  values.  However  at  Separation  and  Vortex-core  stations 
dissipation  is  progressively  increased  compared  to  Station  5.  Pressure-strain  values 
decrease  in  magnitude  progressively  at  Station  5  and  at  Separation  station  compared  to 
2-D  values.  At  the  Vortex-core  station  the  sign  is  changed  to  positive  and  the  magnitude 
is  as  high  as  in  2-D  flow.  Large  positive  peak  value  of  convection  is  shifted  away  from  the 
wall  progressively  at  2-D,  Station  5,  and  Separation  stations.  The  Vortex-core  value  is 
large  very  near  the  wall  but  shows  a  larger  negative  value  away  from  the  wall  at  y'^^SO. 
Pressure  diffusion  is  practically  zero  for  this  stress  only  with  slightly  elevated  values  at 
the  Vortex-core  station.  Viscous  diffusion  term  stays  on  the  same  order  of  magnitude  for 
all  stations.  The  existence  of  a  large  vortex  in  the  outer  region  alters  the  variation  of  the 
near  wall  flow  entirely  at  the  Vortex-core  station.  Pressure-strain  term  becomes  positive. 
Convection  term  is  negative  and  turbulent  diffusion  term  is  as  high  as  the  convection 
term.  Production  is  reduced  and  becomes  negative  around  y*=60.  Dissipation  term 
variation  is  close  to  the  Separation  station  value. 


The  effect  of  three-dimensionality  on  the  'F  transport  is  more  pronounced  then  the  IF 
transport.  Compared  to  2-D  station,  at  the  3-D  stations  both  the  convection  and  turbulent 
diffusion  terms  are  reduced  near  the  wall.  At  the  Vortex-core  station  both  the  convection 
and  turbulent-diffusion  profiles  show  large  negative  peaks  away  from  the  wall  around 
y«200.  Dissipation  peak  value  is  increased  with  increased  three-dimensionality.  The 
pressure-strain  peak  values  near  the  wall  are  increased,  and  the  effect  of  the  term 
extends  further  in  to  the  layer  with  increased  three-dimensionality  and  shows  two  maxinia. 
The  pressure-diffusion  peak  value  located  at  y*=20  for  all  the  cases  increase  with 
increased  3-D  effects.  For  the  Vortex-core  station  the  convection,  production,  and 
turbulent  diffusion  terms  become  important  away  from  the  wall. 


For  a  2-D  TBL  the  production  of  transport  rate  is  zero  which  results  in  the 
pressure-strain  and  the  dissipation  to  be  the  major  terms  in  the  transport.  However  for  the 
3-D  flows  the  production  term  is  one  of  the  major  terms.  The  dissipation  values  gradually 
increase  with  increased  three-dimensionality.  Convection  is  reduced  towards  the  Vortex- 
core  station.  The  term  at  the  Separation  station  is  practically  zero  and  at  the  Vortex-core 
station  the  convection  becomes  negative  with  a  minimum  around  y'^=45.  Turbulent 
diffusion  negative  minima  above  y'^^IO  decrease  with  increased  three-dimensionality. 
Pressure-strain  terms  at  with  positive  values  at  some  regions  decrease  and  become 
negative  through  out  the  layer  towards  the  Vortex-core  station.  The  magnitude  of  the 
pressure-strain  term  is  highly  increased  more  than  two  times  by  the  Vortex-core  station 
compared  to  the  Separation  station. 

Uv 

Three-dimensionality  of  the  flow  results  in  less  in  production  magnitude  and  larger 
in  magnitude  pressure-diffusion  term  at  Station  5  compared  to  2-D  flow.  Towards  the 
Vortex-core  station  the  production,  the  pressure-diffusion,  and  the  pressure-strain  term 
magnitudes  increase.  Even  though  the  variations  are  confined  to  a  narrower  region  near 
the  wall,  the  variations  of  the  profiles  are  similar  to  each  other  in  the  near  wall  region 
including  the  Separation  station.  At  the  Vortex-core  station  the  effect  of  the  outer  region 
vortex  results  in  large  negative  production  and  large  positive  pressure-strain  terms  in  the 
outer  region.  The  y*  location  where  pressure-strain  and  pressure-diffusion  terms  are  equal 
very  near  the  wall  shifts  slightly  away  from  the  wall  at  the  vortex  core  station. 

Uw 

Both  at  the  Separation  and  Vortex-core  stations  the  negative  production  and  the 
positive  pressure-strain  values  are  more  than  twice  of  the  Station  5.  At  the  Vortex-core 
station  this  is  accompanied  by  the  positive  production  and  the  negative  pressure-strain 
values  away  from  the  wall.  Convection  becomes  important  only  at  the  Vortex-core  station 
away  from  the  wall  at  y'=40  to  400  region. 

viv 

Pressure-strain  term  shows  double  peaked  variation  at  the  Separation  station  and 
the  triple-peaked  variation  at  the  Vortex-core  station.  At  the  Separation  and  Vortex-core 
stations  production  profiles  show  two  peaks  with  opposite  signs  to  each  other.  Pressure- 
diffusion  term  for  the  Vortex-core  station  shows  two  peaks  one  at  y'^IO  and  the  other  at 
yWO,  which  corresponds  in  magnitude  to  the  one  observed  at  y^'^IS  to  20  for  the  Station 
5  and  Separation  stations. 


CHAPTER  III  -  THEORETICAL  PRESSURE-STRAIN  MODELS 


In  this  Chapter  seven  theoretical  pressure-strain  models  reported  in  the  literature 
are  briefly  described.  The  models  belong  to  Gibson-Launder  (1978),  Launder-Reece-Rodi 
(1975),  Oberlack-Peters  (1993),  Fu-Launder-Tselepidakis  (2  models),  Shih/Lumley- 
Choi/Lumley,  Speziale-Sarkar-Gatski  (1991). 

Pressure  fluctuations  at  a  point  in  the  flow  can  be  expressed  by  taking  the 
divergence  of  the  Navier-Stokes  equations.  Following  Chou  (1945)  the  equation  for  the 
pressure  fluctuations  is  written  as: 


.  =2  0  U  +11  U  +U  U  ■ 

The  customary  means  of  splitting  this  equation  is  1)  "Rapid"  pressure  term,  which 
is  the  first  term  on  the  right-hand  side  of  the  equation,  and  2)  "Slow"  term  (1978,  Lumley), 
or  the  return  pressure  (Mansour  et  al.,  1988)  which  are  the  last  two  terms.  The  pressure- 
strain  term  obtained  using  the  pressure-fluctuation  equation  is  also  split  into  two  parts  and 
these  parts  are  named  similarly  as  the  "rapid"  term  which  involves  the  mean  velocity 
gradients  and  the  "slow"  term.  The  "slow"  term  is  also  called  the  "return-to-isotropy"  term 
since  it  drives  the  turbulence  into  isotropic  state  (Shih  et  al.,  1993). 

In  the  pressure-strain  model  equations  the  following  abbreviations  are  used. 


6,y=1  when  i=j 
=0  when 
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^  =  U^U^,  £= -  ,  k=—U^U^ 

where  e,y  is  the  dissipation  tensor,  is  the  rate  of  strain  tensor,  £2,^  is  the  vorticity 
tensor,  and  e  is  the  dissipation  of  the  turbulent  kinetic  energy  "TKE=k".  (,i)  denotes  partial 
derivative  of  the  quantity  along  the  "i"  th  axis. 

In  each  model  the  first  line  is  the  model  for  the  "return-to-isotropy"  part,  other  lines 
are  for  the  "rapid"  part  of  the  pressure-strain  (redistribution)  term. 


Ill-A)  Gibson-Launder  model  (1978) 


The  Gibson  and  Launder  (1978)  model  was  used  to  calculate  the  gravitational 
effects  on  the  atmospheric  boundary  layer  and  on  a  nearly  homogenous  free-shear  flow. 
The  model  is  linear  in  Reynolds  stresses  and  similar  to  Launder-Reece-Rodi  model.  This 
model  was  used  by  Schwarz  and  Bradshaw  (1994)  as  Launder-Reece-Rodi  model. 


+C^[(U; 
where,  Ci=3.6,  C2=0.6 


+C2  [  (  Uj  ,  +  UjU,  U.  ,) 


Gibson-Launder  model  has  a  near  wall  correction,  additional  terms  both  for  the  return  to 
isotropy  and  rapid  terms.  The  correction  for  the  return-to-isotropy  was  adopted  from  Shir 
(1973),  and  the  correction  for  the  rapid  term  was  developed  following  this  model  (Gibson, 
Launder,  1978). 


and 


u,Ujn,n^)  f{-—) 
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%,2w  =  ^2 
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where 


Ci'=0.5  ,  Co  =0.3  ,  n.=(  0,1,0)  ,  f=—  ,  l=.LJHL-  ,  k=0.41  and  r  is  the  position  vector. 
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For  consistency  only  one  near-wall  correction  method  was  applied  to  all  the  models.  The 
model  used  is  the  Launder-Reece-Rodi  (1975)  near  wail  correction  model  given  below. 


Ill-B)  Launder-Reece-Rodi  model  (1975) 

One  of  the  earlier  models  presented  for  the  pressure-strain  model  is  of  Launder, 
Reece,  and  Rodi,  (1975).  The  model  was  tested  against  free  shear  flows,  plane  channel 
flow,  and  distorted  homogenous  flows  data  with  reasonable  success.  They  concluded  that 
the  discrepancies  between  computations  and  the  data  were  due  more  to  the  modelling 
of  dissipation  and  diffusion  processes  or  to  the  uncertainties  in  the  measured  quantities 
rather  than  the  pressure-strain  modelling.  Return  to  isotropy  term  was  adopted  from 
Rotta,  (1951).  Rapid  term  away  from  a  boundary  was  developed  using  another  Rotta 

(1951)  approximation  that  ^^pi^={dU,/dxJ  alj’’,  where  a,™  is  the  most  general  form  of 
a  fourth  order  tensor  which  satisfies  a  series  of  kinematic  constraints. 


<I>  =  -c.eb 


+2'^[0.2Sjj+ - ^ - 22 - ^^Ik^jk'^^jk^ik)] 


where,  c^=3.0,  C2=0.4. 

The  model  also  contains  an  additive  near  wall  correction  to  the  pressure-strain  correlation. 


O,, [ 0.125  1  (uiur^kd^j)  -0.015(P,-D,)  ] 
s  y  K  o 

where  y  is  the  distance  perpendicular  to  the  wall.  This  model  was  used  by  Shih  and 
Lumley  (1993)  as  Launder-Reece-Rodi  model. 

Ill-C)  Oberlack-Peters  model  (1993). 

Oberlack  and  Peters  developed  their  pressure-strain  model  based  on  the  two-point 
correlation  function  transport.  This  correlation  function  was  asymptotically  expanded  and 
the  leading  order  terms  were  kept.  This  resulted  in  a  mathematically  simplified 
presentation  of  the  correlation  function.  This  function  was  used  to  calculate  the  rapid  part 
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of  the  pressure-strain  correlation.  The  model  is  used  here  together  with  Rotta’s  (1951) 
retum-to-isotropy  formulation.  The  model  originally  does  not  include  a  near  wall 
correction. 

O  =  -c.eb 

ij  fj 

0.2  s,  .o.6(ft,,s,.,  -|s,,b„s,o  *0.2 

where  Ci=3.0.  This  model  is  similar  to  Launder-Reece-Rodi  model  in  form  however  the 
constants  are  different. 


Ill-D)  Speziale-Sarkar-Gatski  model  (1991) 


SSG  model  originally  does  not  include  a  near  wall  correction  for  the  pressure-strain 
modelling  and  uses  the  homogenous  turbulent  flow  approximations  throughout  the 
mathematical  simplifications.  With  homogenous  flow  assumption  for  turbulent  flows  at  high 
Reynolds  numbers  the  turbulent  diffusion  and  pressure-diffusion  terms  drop  out  and  the 
dissipation  can  be  assumed  isotropic.  The  model  was  developed  to  satisfy  many  physical 
and  asymptotical  consistency  constraints.  Model  results  were  compared  to  LRR  model 
results  for  different  types  of  flows  with  different  numerical  simulations.  Large-eddy 
simulation  and  direct  numerical  simulation  results  show  better  performance  by  the  SSG 
model.  The  model  was  calibrated  using  data  from  isotropic  decay,  return  to  isotropy,  and 
homogenous  shear  flow  experiments  along  with  Rapid  Distortion  Theory  results. 


4.„=-e[c,  6,-3(c, 


c.  o  c* 

^  f  u  o  .  u  c<  ^sko\.  5 


Ij  ^  ^  ^jk ^jk ^ik  -^^ij^kl^kl) 


Where  033=03^, 


( ^ij^ij)  ’P=- ^i^j ,  c,  1  =  1 .8 ,  C31  =  1 .3 ,  =  1 .25 ,  Cg  =0.4 


lll-E)  Shih-Lumley  (1985),  Choi-Lumley  (1984)  model 

*29^[0.2S,*3a,(b,.S,,*i>,,S,.-|8„b„S„) 

+—{2-7  a^)  b^f.Clfi^)  +0.2{b^jjSi,+b^j,Sji-2bi^jb,jSi^i-3bjjbi^iSi^i) 

+o.2{b%a„*b%a„)] 

If  III  >  0  then 
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|3=2.-£ 


p1/2 


1.GX^ 


Q*  p1/2  Q 

y=Jl - 

1+Gx"  ^ 

X=^,  G=-x‘+0.8x° 

Re^ll.  //=-^,  ///=^ 

9£v  2  3 


If  III  <  0  then 

p=2+-^exp[-ZlIZ.][_I£_+80.1  ln(1  +62.4(-//+2.3///))] 


Re 


Y=0 

Re=J‘ 


9ev 


Following  parameters  are  not  dependent  on  the  sign  of 


a. 


=4(1-^0.8F“),  F=1  +9t>„fc,.A,-|‘>A 
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m-F)  Fu-Launder-Tselepidakis  model  I  (1987) 

The  rapid  part  of  the  Fu,  Launder,  Tselepidakis  (FLT)  model  was  developed 
following  Shih  and  Lumley’s  (1985)  approach  Shih  and  Lumley  kept  the  quadratic 
products  of  the  stress  anisotropies.  FLT  kept  including  the  cubic  products.  A  similar 
approach  for  the  return-to-isotropy  term  was  also  developed,  (Launder,  1988).  Shih  and 
Lumley,  (1993  )  point  out  that  this  model  does  not  satisfy  the  Schwarz’s  inequality 
between  the  velocity  and  scalar  fields. 

.8  /I  “  H  £  [b,  ^  1 .2  ((>,„  1>„,  -is,  ll)]-b,it 

-0.6{P,i-h„PJ*0.3tb,i{^) 

-r[-2  ;/(P„-D„) 
where 


?1 


/^i »  n  r%T-Err^  Ml 


k=hi;u. 

2  '  _ 

Dij=-{u;Lr,u,j+UjU,u^^) 

r=0.65,  yA  =  1 -|.(-2//-3///) 

This  model  was  used  as  Fu-Launder-Tselepidakis  model  by  Schwarz  and  Bradshaw 
(1994). 


Ill-G)  Fu-Launder-Tselepidakis  model  II  (1987) 


4.,, =8  “  //  e  [b,  *  1 .2  -  |5,5  II )  ]  -b„e 

+2?'  I  0.2S„t0.3(b„S(,tb„S„-i5,b„S„) 

+0.2 {b^ii  Sji+b^ji  Sji-2 bi^jbijSi^iS bjjbi^iSi^i) 

+0.2{b^i,Qji+b^ji  a,.,) 

+r[4b2„„  [b^i^Q.j„+bji^Q.J+'\2b^jb^j{b^i^Q.„i,+b„^Q.„i,)]  ] 

where 

r=0.7,  >A  =  1-|(-2//-3///) 

This  model  was  used  as  Fu-Launder-Tselepidakis  model  by  Shih  and  Lumley  (1993). 


Ill-H)  Lai-So  model  (1990) 

Lai  and  So  chose  to  model  the  velocity-pressure-gradient  term,  which  includes  both 
the  pressure-strain  and  pressure-diffusion  terms  together.  They  assumed  that  the  existing 
pressure-strain  models  work  satisfactorily  away  from  the  wall  since  they  agree  we  I  with 
the  data  in  that  region.  The  modelling  effort  was  directed  to  formulate  an  additional  terrn 
which  will  not  only  correct  the  pressure-strain  model  for  the  near  wall  effects  but  also  wil 
model  the  pressure-diffusion  term.  The  model  consists  of  an  anisotropic  dissipation  model 
similar  to  that  of  Launder  and  Tselepidakis  (1990)  model  and  a  new  near  wall  c^rrecbon 
similar  to  that  of  Gibson-Launder  (1978)  and  Launder-Reece-Rodi  (1975)  models.  The 
model  adopts  the  Launder-Reece-Rodi  model  as  the  pressure-strain  model  that  needs 

to  be  modified  near  the  wall. 

The  additional  term  to  the  LRR  model  for  the  right  hand  side  of  the  equation 
is  given  by: 

£ 
k 

_ dU:  _ dU, 

where 

uy^k 

to  be  used  together  with 


O,. 


ij,  w 


,  P=-iP„  ,  a’-0.45  ,  n.=(0,1,0) 

3x,,  uXt  2 


C,=1.5 
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^  + u^u,n,n^^n^nju,u,n,n,] 


■  u,u,n,n, 


where  =exp[-( — L)^]  ,  R^= —  The  model  did  not  have  any  observable  effect  on 

150  VE 

the  near  wall  numbers  of  the  2-D  flow,  even  with  different  forms  for  f„  Other  different 
forms  that  were  examined  included  =exp( -f?7./20)  and  f^  ^=exp{-Rj/60),  which 
were  also  tested  by  Launder  and  Tselepidakis  (1990).  The  f,^  ^  varies  between  0.1  to 
0.001  for  different  forms  and  additional  term  effects  make  barely  any  observable  changes 
below  y'^^IO.  The  model  performs  very  similar  to  the  LRR  model.  Therefore  this  model 
is  excluded  from  the  comparisons. 
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CHAPTER  IV  -  COMPARISON  OF  THE  THEORETICAL  AND 
EXPERIMENTAL  PRESSURE-STRAIN  TERMS 


IV-A)  No  near  wall  correction  and  isotropic  dissipation 
approximation 

IV-A-1 )  2-D  flow  station 

Except  the  FLT  models  all  the  models  and  the  data  follow  each  other  closely  above 

y+=40  for  the  TF  transport  pressure-strain  values  (Fig.  57).  For  the  IF  only  GL,  LRR, 
SLCL,  and  SSG  model  predictions  agree  with  the  data  above  y"»40  (Fig.  58).  Near  wall 

variation  of  pressure-strain  term  up  to  y‘^=25  is  well  predicted  by  the  FLT1  model, 

above  y^=35  however  LRR,  GL,  SLCL,  SSG  models  do  a  superior  job  (Fig.  59).Tiv 
pressure-strain  term  estimated  by  the  SLCL  model  above  y‘^«20  is  closest  to  the  data 
values  (Fig.  60).  For  all  the  pressure-strain  terms  the  data  values  and  the  predictions 

above  y"^®!  00  are  close  to  zero  and  close  to  each  other.  Very  near  wall  values  for  IP,  P, 
and  uv  are  not  very  well  predicted  by  the  models,  uw  and  vw  pressure-strain  values  and 
predictions  are  given  in  figures  61  and  62. 

IV-A-2)  Station  5 

The  IP  pressure-strain  term  above  y''=50  is  well  predicted  by  all  the  models  accept 

the  FLT  models  (Fig.  63).  All  the  models  underpredict  the  value  very  near  the  wall.T^ 
pressure-strain  term  is  very  well  predicted  by  the  OP ,  SLCL,  and  FLT1  models  till  y  ~8 
(Fig.  64).  Between  y‘^=8  to  60  none  of  the  models  follow  the  trend  of  the  data  values. 
Above  y*=60  GL,  LRR,  SLCL,  and  SSG  values  follow  the  data  values  closely.  The 

variation  of  PP  is  not  well  predicted  throughout  the  layer  except  in  y^»30  to  60  region  by 
FLT2  model  (Fig.  65).  Although  the  trend  of  the  predictions  and  the  trend  of  the  data  are 

same  both  values  are  close  to  zero.  77v  pressure-strain  values  above  y'^^SO  are  equally 
well  predicted  by  the  LRR,  GL,  and  SSG  models  (Fig.  66).  Above  y^=60  all  the  model 
estimates  are  close  to  the  data  and  close  to  the  zero  value.  Near  the  wall  between  y^«1 5 
to  30  SLCL  model  values  are  very  close  to  data  values.  Trend  of  the  data  nearest  the 

wall  is  not  followed  by  any  of  the  models,  uiv  pressure-strain  term  was  well  predicted 
above  y^=15  by  the  FLT1  model  (Fig.  67).  Near  wall  prediction  of  the  model  is  also  not 
far  away  from  the  data.  OP,  SSG,  LRR,  GL  models  do  a  better  job  below  yW.Between 

y*=15  to  30  only  SLCL  model  predicts  the  Pw  pressure-strain  variation  (Fig.  68).  Above 
y^=40  all  model  estimates  and  the  data  values  are  both  close  to  each  other  and  to  the 
zero  value. 


IV-A-3)  Separation  station 


TP  pressure-strain  value  is  not  well  estimated  by  any  model  at  this  station  except 

when  the  data  values  and  estimates  are  very  close  to  zero  (Fig.  69).  pressure-strain 
variation  is  well  approximated  only  above  y*=40  (Fig.  70).  Even  though  the  estimates  and 

the  data  values  overlap  at  certain  regions  none  of  the  models  were  able  to  predict  the 
data  values  until  y"’~500  (Fig.  71).  uv  value  is  predicted  well  by  the  SLCL  model  above 

y‘'=20  (Fig.  72).  The  trend  of  SLCL  and  FLT1  models  are  closer  to  the  trend  of  XheUw 
data  throughout  the  layer,  with  superior  agreement  of  FLT1  model  (Fig.  73).  Above  y'^^SO 

all  model  values  and  data  values  are  close  to  each  other,  vw  pressure-strain  values  are 
equally  well  predicted  by  the  FLT1  and  SLCL  models  above  y'^IS  (Fig.  74). 

IV-A-4)  Vortex-core  station 


The  variations  of  TP,  T^.  W  are  not  well  predicted  by  any  model  (Figures  75,  76, 
and  77).  Values  predicted  by  FLT  models  above  y‘'=300  are  closer  to  the  pressure-strain 
data  values.  Shear  stress  pressure-strain  estimates  at  this  station  are  also  much  different 
than  the  data  values  through  out  the  layer  especially  near  the  wall.  Estimated  values 

follow  the  same  trend  of  the  data  above  y'-l  00,  y"’=400,  and  y‘"=200  for  the  Uv,  TJw,  and  vw 
stresses  respectively  (Figures  78,  79,  and  80). 

IV-B)  Algebraic  anisotropic  dissipation  rate  models 

The  fact  that  the  limiting  value  of  the  dissipation  at  the  wall  is  not  isotropic  is 
realized  by  many  researchers  and  different  algebraic  models  are  proposed  in  the 
literature.  Nagano  and  Tagawa  (1991)  argued  that  simple  algebraic  models  could  not 
capture  the  physics  and  developed  a  new  closure  scheme  which  is  not  algebraic.  In  this 
report  several  algebraic  dissipation  rate  models  are  summarized.  For  consistency  only 
one  model  (Hallback,  et  al.,  1990)  was  used  and  tested  in  comparison  of  measured  and 
predicted  pressure-strain  models. 


IV-B-1)  Launder  and  Reynolds  (1983)  model 


Launder  and  Reynolds  (1983),  using  the  near  wall  variations  of  the  fluctuating 


U:U, 


velocity  components,  show  that  the  common  assumption  of  e.^.=(_Li)e  underestimates 

k 


the  dissipation  in  the  v^  transport  equation  by  a  factor  of  4.  They  propose  a  model  which 
satisfy  the  near  wall  variation: 


where  nj=(0,1,0),  a  =  1,  and  P=1  with  the  provisions  that  for  or  £33,  TP  is  negligible 

compared  to  TP  or  Tip. (Therefore,  |3=0  when  i=1,j=1  and  i=3,j=3).  This  equation  shows 
that  the  dissipation  of  each  shear  stresses  is  different  than  zero. 


IV-B-2)  Hallback,  Groth,  and  Johansson  model  (1990) 
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Hallback,  Groth,  and  Johansson  developed  an  explicit  model  relating  the 
anisotropies  of  the  dissipation  rate  tensor  to  the  anisotropies  of  the  Reynolds  stress 
tensor  using  Lumley’s  invariant  theory  (1978).  The  general  expression  for  the  dissipation 
rate  anisotropies  given  by  the  invariant  theory  was  expanded  in  power  series  of  the 
Reynolds  stress  anisotropies.  A  series  of  mathematical  and  physical  constraints  were 
employed  on  this  expansion  to  obtain  the  dissipation  rate  anisotropy  tensor  ey. 

-j)  -|8,/  a"'* 

~uir,  2  c  ,1 _  ..  3 

where,  a,y=(_^) --5,^  ,  l<=-u.,Ui  ,  //=a,.^a^, ,  a=- 

The  model  was  tested  against  numerical  simulations  of  homogeneous  turbulence 
with  success. 

IV-B-3)  Lai,  and  So  model  (1990) 

Lai  and  So  developed  an  anisotropic  dissipation  rate  tensor  model  similar  to  the 
Launder  and  Reynolds  (1983)  model.  The  model  does  not  need  any  provisions  for 
and  833. 


where 


C.1 


1  +^u^u,n,uj{2k) 


C.i=exp[-(R^/150)2]  ^ 
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IV-B-4)  Launder,  and  Tselepidakis  model  (1990) 


8;^=4e*+(1 

where 

'  ^  5  , 

r  1  4-  ^  p  q  p  q-l 

^  2  k 

and  f  =exp[-R/20]  ,  Rf=—  ,  k=^lJ^j 

IV-C)  Anisotropic  dissipation  rate  effects  on  the  experimental 
stress  transport  rate  budget 

IV-C-1)  Anisotropic  dissipation  rate  distribution 
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Anisotropy  of  the  dissipation  rate  tensor  affects  the  pressure-strain  rate  term 
extracted  from  the  Reynolds’  stress  transport  equations.  Calculation  of  convection, 
production,  and  viscous  diffusion  are  not  affected  by  the  anisotropy  of  the  dissipation  rate 
tensor.  Since  the  pressure-strain  term  is  extracted  from  the  balance  of  the  budget  of 
transport  equations,  it  is  the  only  term  affected.  Anisotropy  of  the  dissipation  term  also 
does  not  affect  the  TKE  transport  budget.  Figures  81  to  84  show  the  distribution  of  the 
dissipation  rates  for  different  stresses.  The  figures  also  include  the  isotropic  dissipation 

rate  for  comparison.  For  the  2-D  case,  figure  (81)  shows  that  while  the  dissipation 

rate  is  very  similar  to  the  isotropic  dissipation  rate,  the  IF  dissipation  rate  is  much  larger 

than  both  isotropic  and  IF  dissipation  rates  (Fig.  81).  While  the  peak  value  of  theu^ 

dissipation  is  close  to  2  times  of  the  isotropic  dissipation  rate,  the  dissipation  in  the 

transport  is  halved.  The  dissipation  in  the  Fw  and  vw  transport  are  practically  zero  and 

the  dissipation  in  the  transport  of  Fv  is  negative  with  a  peak  value  of  approximately  1/3 

of  the  'F  dissipation. 

At  Station  5  the  anisotropy  of  the  TF  dissipation  rate  is  noticeable.  Compared  to 

the  2-D  case  the  F  dissipation  rate  is  reduced  (Fig.  82).  The  Fw  dissipation  rate  is 

higher  in  magnitude  than  the  F  dissipation  rate.  The  Fv  and  Fw  dissipation  rate 

absolute  values  reach  to  half  of  the  F  dissipation  rate,  however  they  are  an  order  of 

magnitude  smaller  than  the  F  and  FF  dissipation  rate  values.  Dissipation  rate  is  close 
to  isotropic  distribution  above  y''«100. 

The  magnitude  of  the  dissipation  rate  values  in  general  are  increased  more  than 

twice  by  the  Separation  station  compared  to  the  2-D  station  values  (Fig.  83).  TheF 

dissipation  rate  is  the  highest  near  the  wall.  Above  y"‘»25,  F  and  FF  dissipation  rate 

values  follow  each  other  closely.  The  F  dissipation  rate  is  relatively  smaller  than 
dissipation  rate  values  compared  to  Station  5.  Shear  stress  dissipation  rate  terms  are  also 
increased. 

At  the  Vortex-core  station  above  y'^=20,  FF  dissipation  rate  is  larger  is  larger  than 
dissipation  rate  (Fig.  84).  The  magnitudes  of  the  dissipation  rate  values  are  higher  than 

compared  to  other  stations  and  difference  between  the  isotropic  distribution  and  F  and Tv^ 

dissipation  rate  values  are  further  increased.  The  Fw  dissipation  rate  changes  sign 
around  y'’'=20  and  shows  a  positive  peak  around  y'’’=40.  Normal  stress  dissipation  rate 
terms  follow  each  other  closely  above  y'=125  and  show  a  relative  maxima  at  y"=300. 
Stress  dissipation  rate  values  above  y^=1 00  are  close  to  zero. 

IV-C-2)  Effects  of  anisotropic  dissipation  rate  on  the 
experimentally  extracted  pressure-strain  terms 

Following  the  Hallback  et  al.  model,  the  dissipation  process  near  a  wall  is  highly 

anisotropic  with  most  of  the  dissipation  occurring  in  the  F  and  the  FF  transport. 
However  above  y'^^lOO  the  dissipation  rate  is  isotropic  for  the  3-D  flows.  In  the  2-D  flow, 
the  dissipation  rate  becomes  isotropic  at  a  higher  y*  (y'’’=350)  location.  The  dissipation 
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rate  values  for  the  IF  and  'W  transport  are  increased  by  the  three-dimensionality  of  the 
flow  and  the  shear  stress  transport  dissipation  rate  values  are  away  from  zero  with  values 

an  order  of  magnitude  smaller  than  the  IF  values. 

The  anisotropy  of  the  dissipation  rate  effects  only  the  pressure-strain  value 
extracted  from  the  balance  of  the  budget  of  a  stress.  The  increase  in  the  dissipation  rate 
value  of  a  stress  above  the  isotropic  distribution  value  at  a  point  in  the  profile  also  results 
to  an  increase  by  the  same  amount  to  the  pressure-strain  value  at  the  same  point  in  the 
profile  extracted  using  isotropic  dissipation  rate  distribution.  This  results  in  reduced  values 

in  magnitude  of  the  pressure-strain  values  for  the  IF  in  the  2-D  flow  (Fig.  85). IF 
transport  pressure-strain  term  becomes  negative  very  near  the  wall  below  y*=20  (Fig.  86). 

The  magnitude  of  the  dissipation  rate  for  the  W  stress  is  too  small  to  have  any 

observable  changes  in  the  pressure-strain  term  of  this  stress  (Fig.  87).  TiF  pressure-strain 
distribution  is  very  similar  to  the  one  calculated  with  isotropic  dissipation  rate 
approximation  (Fig.  88). 

At  Station  5  location  the  effect  of  anisotropic  dissipation  rate  in  the  IF  transport 

is  to  reduce  the  magnitude  of  IF  pressure-strain  rate  (Fig.  89).  The  values  near  the  wall 
below  y*~^  5  become  close  to  zero.  Negative  peak  value  at  y'"=20  is  halved.  The  pressure- 

strain  profile  of  IF  transport  shows  two  distinct  positive  peaks,  one  at  y'"=8  and  one  at 
y^=50,  and  a  negative  peak  at  y^=20  (Fig.  90).  The  peak  dissipation  rate  value  is  2/3  of 

the  positive  peak  value  of  the  pressure-diffusion  term.  Effect  on  liF  pressure-strain  term 
is  felt  at  y*~25,  the  positive  peak  value  is  increased  by  0.01  from  a  value  of  0.015  (Fig. 

91).  The  Uv  and  vw  pressure-strain  distributions  (Figures  92,  and  93)  are  not  effected 
since  the  dissipation  rate  predicted  for  these  stresses  are  small  in  magnitude.  However 

the  peak  value  of  the  pressure-strain  term  for  the  Uw  transport  is  decreased  by  -0.035 
from  a  value  of  -0.125  (Fig.  94). 

At  the  Separation  station  the  effect  of  the  anisotropy  results  in  positive  pressure- 

strain  values  near  the  wall  below  ^=20  for  the  IF  transport  equation  (Fig.  95).  The  profile 
shows  a  positive  peak  at  ^=20  rather  than  a  negative  peak  at  y*«10  that  was  observed 
with  isotropic  dissipation  rate  approximation.  The  profile  for  the  pressure-strain  term  of 

the  IF  transport  is  also  drastically  affected  (Fig.  96).  This  profile  shows  a  negative  peak 
at  y'’=18  rather  than  a  positive  peak  around  y^=10  with  reduced  values  throughout  the 
layer.  The  variation  becomes  a  mirror  image  of  the  pressure-diffusion  around  a  value  of 

0.01 .  The  values  of  the  pressure-strain  in  the  transport  equation  of  UF  (Fig.  97)  are  about 
constant  between  y*«1 0  to  80.  Both  the  Uv  and  vw  pressure-strain  term  profiles  are  very 

similar  to  the  profiles  obtained  with  isotropic  dissipation  rate  (Figures  98,  and  99).Uw 
pressure-strain  profile  is  also  similar  with  slightly  reduced  peak  value  (Fig.  100). 

At  the  Vortex-core  station  the  pressure-strain  values  of  the  IF  transport  are 
elevated  below  y^=100,  resulting  in  larger  values  than  the  production  values  in  y^«15  to 

30  range  (Fig.101).  The  variation  of  the  IF  pressure-strain  rate  is  also  highly  affected  by 
the  anisotropy  of  dissipation  rate  as  in  the  other  stations  (Fig.102).  While  the  distribution 
is  positive  till  y*«400  where  the  dissipation  rate  is  isotropic,  the  profile  shows  a  negative 
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peak  at  y*~2S  with  anisotropic  dissipation  rate  approximation.  The  pressure-diffusion  term 

becomes  one  of  the  most  important  terms  below  y‘'=50  in  describing  the  transport  of  17^. 

The  pressure-strain  term  is  increased  by  an  amount  of  0.1  below  y"=50  (Fig. 103).  The 
variation  is  similar  to  results  when  the  isotropic  dissipation  rate  is  used.  Pressure-strain 
terms  in  the  shear  stress  transport  equations  are  not  affected  by  the  anisotropy  of  the 
dissipation  rate  since  the  anisotropic  dissipation  rate  values  for  this  station  are  close  to 
zero  (Figures  104,  105,  and  106). 

IV-D)  Anisotropic  dissipation  rate  effects  on  theoretical  and 
experimental  pressure-strain  comparison 

Anisotropic  dissipation  rate  approximation  of  Hallback  et  al  (1 990)  does  not  affect 
the  pressure-strain  distribution  estimated  by  the  models.  However,  it  affects  the 
experimental  values  extracted  from  the  stress  transport  budget  balances.  Anisotropic 
distribution  results  in  better  agreement  of  the  predictions  and  the  data  of  different  models 
at  different  regions.  In  the  next  paragraphs,  anisotropic  dissipation  rate  effects  in 
predictions  of  the  models  will  be  discussed. 

IV-D-1)  2-D  Station 


The  IF  pressure-strain  term  is  well  predicted  by  the  FLT2  model  near  the  wall 
below  y'=25  (Fig.  107).  None  of  the  models  predicts  the  data  in  this  range  with  the 
isotropic  assumption.  Between  y"^=25  to  50,  none  of  the  models  predicts  the  variation; 
above  y'"=50  all  the  model  estimates  follow  each  other  and  the  data  values  closely  except 

the  FLT  models.  The  IF  pressure-strain  is  also  well  predicted  by  the  FLT2  model  near 
the  wall  below  y''=25  (Fig. 108).  Above  y''==40,  data  values  are  closer  to  FLT1  and  OP 

model  predictions  rather  than  the  GL,  LRR,  SLCL,  or  SSG  models.  Both  the  TF  and  the  tTv 
pressure-strain  values  are  similar  to  the  isotropic  dissipation  rate  approximation  variation 

since  the  dissipation  rate  distributions  for  TP  and  Uv  are  very  similar  to  that  of  isotropic 
dissipation  rate  (Figures  109,  and  110). 

IV-D-2)  Station  5 

Pressure-strain  data  near  the  wall  for  the  TF  transport  equation  is  close  to  zero  till 
y^'^l  5.  At  y'^^l  5  the  estimated  profiles  show  negative  peaks  (Fig.1 11).  Above  y''=50  all  the 
model  values  are  not  far  away  from  the  data.  Below  y'=40,  data  and  the  predictions  are 

further  apart  from  each  other  then  the  isotropic  dissipation  rate  approximation.  The'P 
data  are  closely  predicted  by  the  FLT2  model  throughout  the  layer  (Fig.1 12).  This  model 
behaved  as  the  worst  model  with  the  isotropic  dissipation  rate  approximation.  Between 

)/*~25  to  y"=60,  the  TP  data  variation  is  closely  followed  by  the  OP  model  estimates 
(Fig.1 13)  rather  than  the  FLT2  model  with  the  isotropic  dissipation  rate  approximation. 
Above  this  location  data  values  and  all  the  predictions  are  close  to  each  other  and  to  the 

zero  value.  The  variations  of  the  pressure-strain  data  for  Uv  and  vw  are  not  different  than 
the  isotropic  dissipation  rate  distribution  case  (Figures  114,  115).  Predictions  by  the 
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models  for  the  Uw  stress  are  less  in  agreement  with  the  data  in  the  y‘^=10  to  30  range 
(Fig.  1 16).  Variations  near  the  wall  are  well  approximated  by  the  SLCL,  FLT1  models.  Also 
above  y"^*50  all  the  model  estimates  follow  the  data  closely. 

IV-D-3)  Separation  station 

The  anisotropic  dissipation  rate  results  in  much  different  pressure-strain  distribution 

in  the  IF  transport  equation  than  the  isotropic  dissipation  rate  approximation  (Fig. 117). 
The  model  variations  are  also  very  much  different  than  data.  While  the  data  show  a 
positive  peak  at  y=20,  the  models  predict  much  larger  negative  peak  at  y*=15.  The  trend 

of  the  data  is  opposite  to  the  trend  of  the  predictions.  The  "iF  prediction  by  FLT2  is  very 
good  throughout  the  layer  (Fig.118),  while  the  other  models  show  a  positive  peak  value 
at  the  location  where  the  data  show  a  negative  peak  value.  The  FLT2  model  with  the 
isotropic  dissipation  rate  approximation  was  the  only  model  which  predicted  values  with 

an  opposite  trend  to  the  data.  The  W  pressure-strain  data  are  not  predicted  well  by  the 
models  (Fig.  1 19).  Although  the  estimates  are  in  the  vicinity  of  the  data,  the  variations  are 

much  different.  The  Uv  data  values  are  closely  predicted  by  the  SLCL  and  FLT1  models 

above  y^=15  (Fig. 120)  similar  to  with  isotropic  dissipation  rate.  The  Uw  pressure-strain 
data  agreement  with  the  FLT1  model  with  an  isotropic  dissipation  rate  approximation  is 
poorer  (Fig. 121).  However  the  model  distributions  are  close  to  that  of  the  data  above 

y'==25.  The  vw  pressure-strain  distribution  is  not  different  than  isotropic  dissipation  rate 
distribution  and  is  well  predicted  by  the  SLCL  and  FLT1  models  above  y'^^IS  (Fig. 122). 

IV-D-4)  Vortex-core  station 

As  previously  observed  for  all  the  normal  and  shear  stresses  the  pressure-strain 
data  are  very  poorly  predicted  by  all  the  models.  The  differences  between  the  data  and 
the  predictions  increase  (Figures  123  through  128). 

IV-E)  Effect  of  the  near-wall  correction  on  the  theoretical 
pressure-strain  models 

Near  wall  corrections  were  calculated  using  the  Launder-Reece-Rodi  (1975)  model. 
For  consistency  the  model  corrections  were  applied  to  the  other  models.  The  near  wall 
corrections  are  large  below  y^«1 00  for  the  2-D,  Station  5  and  Separation  stations  (Figures 
129,  130,  and  131).  The  corrections  span  up  to  y'"=300  at  the  Vortex-core  station 
(Fig.  132).  The  near  wall  corrections  affect  only  the  mathematical  pressure-strain  models. 
The  corrections  are  additive  to  the  model  estimates.  At  all  3-D  stations  the  largest 

corrections  are  estimated  for  the  Uv  and  Vw  stresses.  At  2-D  station  the  Uv  correction 

is  the  largest.  The  IF  corrections  at  all  stations  are  close  to  one  another.  TheTF 

corrections  are  on  the  same  order  of  IF  corrections  and  become  large  at  the  Vortex-core 

station.  The  IF  correction  is  small  for  the  2-D,  Station  5  and  Separation  stations  and  is 

on  the  same  order  of  magnitude  of  the  IF  correction  for  the  Vortex-core  station.  TheUw 


correction  is  on  the  same  order  of  the  correction  for  the  Station  5  and  Separation 
stations.  However,  at  the  Vortex-core  station  the  correction  is  larger. 


IV-F)  Evaluation  of  near-wail  correction  effects  on  theoretical 
and  experimental  pressure-strain  with  isotropic  dissipation  rate 
approximation 

IV-F-1 )  2-D  Station 

At  the  2-D  station  near-wall  corrections  to  the  normal  stress  pressure-strain  models 
improve  the  prediction  of  the  data  values.  However  differences  between  the  data  and  the 
predictions  appear  above  y'’=40.  Below  ^=40,  the  OP  model  estimates  are  closer  to  the 

data  for  IF  (Fig. 133).  For  IF  and  IF  LRR  and  GL  models  are  relatively  closer  to  the 

data,  respectively  (Figures  134,  and  135).  For  the  Uv  stress  the  corrections  near  the  wall 
below  y'^=20  result  in  highly  overpredicted  values  (Fig. 136). 

IV-F-2)  Station  5 

Near-wall  corrections  affect  the  estimated  values  by  the  models  to  be  closer  to  the 

data  values  for  the  pressure-strain  term  of  IF  (Fig.  137).  However  the  estimates  are  still 
very  much  different  than  the  data.  The  FLT1  model  trend  is  closest  to  the  data  variation. 

Double-peaked  pressure-strain  term  variation  for  the  IF  is  not  followed  by  any  of  the 
models  (Fig.  138).  Near  wall  correction  results  in  reduced  agreement  throughout  the 

profile,  especially  near  the  wall.  Model  predictions  for  the  pressure-strain  variation  of  mF 
are  closer  to  the  data  values  (Fig. 139),  but  the  change  in  the  predicted  values  compared 

to  the  no-near-wall-corrected  estimates  is  very  small.  Near  wall  corrections  for  them/ 
shear  stress  pressure-strain  term  are  opposite  to  the  trend  of  the  data  near  the  wall 
especially  below  y^'^IS  (Fig. 140).  The  pressure-strain  term  is  well  predicted  both  by  the 
GL  and  LRR  models  above  this  \/*  location.  The  agreement  between  the  data  and  all  the 

models’  estimates  except  the  FLT  models  is  better  above  y’'»30forthe  Uw  (Fig.141).  Them 
pressure-strain  term  is  well  predicted  by  the  FLT1,  LRR,  and  GL  models  above  y^=20 
(Fig.142). 

IV-F-3)  Separation  station 

Near-wall  corrections  affect  the  model  estimates  in  a  similar  way  as  at  Station  5. 

The  correction  causes  the  trend  of  the  predictions  for  the  IF  and  IF  normal  stresses  to 
be  closer  to  the  data  variations,  but  the  amount  of  correction  does  not  result  in  collapse 
of  the  data  and  predicted  profiles  (Figures  143,  and  144).  The  near-wall  corrected 

estimates  for  the  HF  results  in  reduced  agreement  (Fig.  145).  The  shear  stress  pressure- 
strain  data  and  predictions  differ  highly  below  y^'^l  5  to  20.  The  Uv  data  profile  is  closely 
followed  by  the  GL,  LRR,  SSG,  and  SLCL  models  above  y''=15  (Fig. 146).  The  Uw  data 
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are  well  predicted  by  all  the  models,  except  the  FLT  models  above  y^«20  (Fig.147).  The  vw 
pressure-strain  estimates  by  the  SSG,  LRR,  and  GL  models  are  in  good  agreement  with 
the  data  above  (Fig. 148). 

lV-F-4)  Vortex-core  station 

The  near-wall  corrections  results  in  good  agreement  of  OP,  FLT2,  and  FLT1 

models  with  the  pressure-strain  term  of  IP  normal  stress  below  y^»15  (Fig. 149).  Data  and 
predictions  of  GL,  SSG,  and  FLT1  are  in  agreement  in  the  yWO  to  150  range.  Also  in 
the  y*=300  to  700  range  FLT  models  predict  the  data  well.  The  data  values  for  the 

pressure-strain  of  do  not  change  in  the  profile  up  to  y*~200  where  the  predictions  in 
the  same  range  highly  vary  (Fig.  150).  In  similar  fashion  as  to  the  values  predicted  without 
the  wall  corrections,  the  predicted  values  with  the  near  wall  corrections  are  far  away  from 
the  data  values,  except  in  the  y*=250  to  450  limited  range.  In  this  range,  the  FLT  models’ 

estimates  are  close  to  the  data  values.  Closer  agreement  of  the  TF  data  and  the 
predictions  in  the  region  above  y^“70  without  the  near-wall  corrections  is  delayed  till 
y"^=150  (Fig.  151).  Above  y‘"=150,  data  and  predictions  of  FLT  models  are  in  good 

agreement.  Near  wall  corrections  bring  the  77v  data  estimates  of  the  GL  model  above 
y‘"»20  and  SLCL  model  above  y'^^SO  closer  to  the  data  (Fig. 152).  The  uw  pressure-strain 

term  is  poorly  estimated  throughout  the  layer  (Fig. 153).  The  vw  data  estimates  of  the 
models  are  opposite  to  the  near  wall  data  (Fig.  154).  Above  y’'«100  the  data  and  the 
predictions  follow  each  other. 

IV-G)  Combined  effect  of  the  near-wall  correction  and  the 
anisotropic  dissipation  rate  on  the  theoretical  pressure-strain  model 
predictions 

It  was  previously  discussed  that  the  anisotropic  dissipation  rate  assumption  results 
in  different  distributions  of  the  experimental  pressure-strain  term  extracted  from  the  stress 
transport  budget.  A  near-wall  correction  term  is  used  to  take  into  account  the  existence 
of  the  wall  in  the  wall  bounded  flows  and  affects  the  pressure-strain  term  calculated  by 
the  model.  In  this  section  the  combined  effects  of  this  correction  and  the  anisotropic 
dissipation  rate  approximation  on  the  agreement  of  the  data  and  the  calculations  will  be 
discussed. 

IV-G-1)  2-D  station 


Even  though  the  IP  pressure-strain  term  is  well  predicted  near  the  wall  below 
y‘’’=20  with  only  an  anisotropic  dissipation  rate  approximation  and  reasonable  well 
estimated  above  y^=15  with  only  the  near-wall  corrections  applied  models,  both  effects 
together  do  not  result  in  better  agreement  (Fig.155).  The  data  are  well  predicted  by  all 

the  models  except  the  FLT  models  above  y^=40.  The  P  pressure-strain  data  were  well 
predicted  by  the  FLT2  using  anisotropic  dissipation  rate  approximation  very  near  the  wall 
(Fig. 156).  The  combined  effect  reduces  this  agreement.  The  region  above  y"'=40  is  well 


estimated  by  different  models  with  or  without  the  correction  or  approximations  made.  The 
pressure-strain  variation  is  very  well  approximated  by  the  models  with  or  without  the 

assumptions  (Fig.  157).  The  Uv  pressure-strain  predictions  are  very  similar  to  the  near¬ 
wall  corrected  predictions.  SSG,  GL,  SLCL  model  estimates  agree  with  the  data  above 
y^*18  (Fig.158). 

IV-G-2)  Station  5 

Estimates  to  the  IF  pressure-strain  term  variation  are  improved  by  the  combined 
effect  below  y*«8  (Fig.159).  However  the  differences  between  the  data  and  the 

predictions  in  the  y"=10  to  50  region  are  not  reduced.  The  combined  effect  on  theT^ 
pressure-strain  term  is  similar  to  the  near-wall  correction  effect.  The  data  and  the 

estimates  differ  highly  below  y''=50  (Fig.160).  The  near  wall  variation  of  IF  data  up  to 
y'"=25  is  not  well  predicted  with  or  without  the  effects  of  dissipation  rate  anisotropy  and 
near-wall  correction  (Fig.  161).  Above  this  location  OP  model  values  follow  the  data 

closely  with  the  current  assumptions.  The  Uv  pressure-strain  variation  is  well  estimated 
by  the  LRR,  GL,  and  FLT1  models  above  y'’=15,  similar  to  near-wall  corrected  model 

estimates  (Fig. 162).  Agreement  between  the  Uw  data  and  the  predictions  appears  to  be 
degraded  with  the  near-wall  correction  or  the  anisotropic  dissipation  rate  approximation, 

or  the  combined  effects  of  these  assumptions  (Fig.  163).  Above  y*=15  the  vw  data  and 
the  estimates  of  LRR,  GL,  and  FLT2  model  agree  closely  (Fig.164).  The  near  wall 
variation  is  not  well  estimated  by  any  assumptions, 

IV-G-3)  Separation  station 


The  variation  of  the  model  predictions  for  IF  pressure-strain  and  the  data  highly 
differ  from  each  other  with  opposite  signed  peaks  at  different  y*  locations  (Fig. 165).  The 

good  agreement  between  the  F  pressure-strain  data  and  the  LRR  model  with  the 
anisotropic  dissipation  rate  approximation  is  canceled  with  the  combined  effect  (Fig. 166). 

The  IF  data  are  well  approximated  by  the  OP  model  in  the  y*^^  8  to  70  region  when  the 

effects  are  combined  (Fig.167).  Combined  effects  results  in  Fv  pressure-strain  data  to  be 
followed  by  the  LRR,  and  GL  models  starting  at  y^=15  (Fig.  168).  Combined  effects  or 

each  effect  singly  degrades  the  agreement  between  the  models  and  the  data  forlTw 

(Fig.169).  For  the  vw  pressure-strain  term  SLCL,  LRR,  and  GL  model  predictions  are  in 
good  agreement  above  y'^=15  (Fig. 170).  However,  the  differences  between  the  variations 
of  the  estimates  and  the  data  are  increased  closer  to  the  wall. 

IV-G-4)  Vortex-core  station 

The  close  agreement  between  the  model  estimates  and  the  data  for  IF  at  this 
station  is  reduced  compared  to  only  near-wall  corrected  model  estimates  (Fig. 171).  The 
differences  between  the  models  and  the  data  for  the  other  normal  and  shear  stress 
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pressure-strain  data  are  very  large  and  adjustments  due  to  the  anisotropic  dissipation  rate 
approximation  and  near-wall  correction  separately  or  together  are  not  large  enough  to 
compensate  for  them  (Figures  172  through  176). 


CHAPTER  V  -  DIRECT  NUMERICAL  SIMULATION  (DNS) 
COMPARISON 

In  this  Chapter  comparison  of  several  existing  DNS  calculations  to  the  present 
transport  rate  budgets  are  made.  DNS  simulations  are  usually  calculated  at  low  Reynolds 
numbers  due  to  very  costly  computer  time  requirements.  DNS  studies  are  also  directed 
to  much  simpler  flow  configurations  than  the  present  study  configuration.  Comparisons 
are  therefore  qualitative  rather  than  quantitative.  General  trends  of  the  transport  budget 
terms  are  outlined  and  compared. 

V-A)  Study  of  Sumitani-Kasagi  (1995) 

Sumitani  and  Kasagi  (1 995)  studied  a  fully  turbulent  channel  flow  with  heat  transfer 
and  uniform  wall  injection  at  a  low  Reynolds  number  (Re=150=channel  half  width  x  / 

v).  Their  results  showed  that  while  the  normal  stresses  and  the  shear  stress  increased 
on  the  injection  side  and  decreased  on  the  suction  side,  skin  friction  increased  on  the 
suction  side  and  decreased  on  the  injection  side.  Balance  of  the  TKE,  normal  stress, 
shear  stress,  and  dissipation  rate  transport  equations  were  calculated  and  presented. 
^The  present  2-D  case  data  of  the  terms  in  the  budget  of  the  transport  equations  show 
similar  characteristics  to  the  presented  DNS  calculations  by  Sumitani  and  Kasagi. 

For  TP  the  similarities  can  be  summarized.  The  production  terms  both  for  the  DNS 
and  the  2-D  data  peak  at  the  similar  y*  locations  on  the  suction  side  (y'^^IS).  Convection 
terms  on  the  suction  side  show  increased  values  below  y"«10.  Pressure-diffusion  is  not 
presented  by  DNS.  For  the  present  data  it  is  calculated  as  approximately  zero.  Pressure- 
strain  terms  on  the  suction  side  are  close  to  the  dissipation  rate  variation  above  y^’^SO 
and  values  are  also  not  far  away  from  the  production  values  with  an  opposite  sign.  A 
similar  variation  is  observed  for  the  present  data.  Dissipation  rate  values  were  observed 
to  be  larger  than  the  pressure-strain  terms  for  the  DNS  study  on  the  suction  side  below 

which  is  more  closely  approximated  by  the  anisotropic  dissipation  rate 
approximation  applied  to  the  present  data.  On  the  suction  side  DNS  results  show  that 
pressure-strain  values  are  close  to  zero,  which  is  not  the  case  observed  with  the  present 
data.  DNS  shows  that  the  dissipation  rate  and  the  viscous  diffusion  terms  are  not  zero 
at  the  wall  for  this  stress. 

Experimental  data  for  Station  5  and  Separation  stations  are  qualitatively  similar  to  the  2-D 
data,  however  the  distributions  for  the  Vortex-core  station  are  very  much  different  both 
from  the  2-D  data  and  the  DNS  solutions. 

The  DNS  results  shows  that  all  the  IP  transport  terms  at  the  wall  are  zero.  The 
pressure-diffusion  term  on  the  suction  side  was  elevated  as  high  as  the  dissipation  rate 
from  near  the  wall  to  about  y‘'«60.  Pressure-strain,  and  pressure-diffusion  terms  on  the 
injection  side  are  equally  high  with  opposite  signs  very  near  the  wall  below  y^'^S, 
convection  was  practically  zero,  viscous  diffusion  was  elevated  near  the  wall  below  y'^IS, 
and  turbulent  diffusion  changed  sign  around  y'^^IS.  Pressure-strain  and  the  dissipation 
rate  terms  are  the  dominant  terms  away  from  the  wall.  Similar  trends  are  observed  for  the 
present  2-D,  Station  5,  and  Separation  stations. 

DNS  results  for  TP  showed  that  the  viscous-diffusion  and  dissipation  rate  terms 
at  the  wall  are  not  zero  for  this  normal  stress.  The  dominant  terms  very  near  the  wall 
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below  ^=5  are  the  viscous  diffusion  and  the  dissipation  rate  terms  both  on  the  suction 
and  injection  sides.  Away  from  the  wall  pressure-strain  and  the  dissipation  rate  term 
values  are  approximately  equal  to  each  other  with  opposite  signs  to  each  other. 
Production  and  pressure-diffusion  terms  were  not  presented  since  they  are  practically 
zero.  These  observations  are  also  valid  for  the  present  2-D  data.  At  3-D  flow  stations  the 
production  term  for  this  stress  is  the  largest  term  and  highly  alters  the  budget. 

Pressure-diffusion  and  pressure-strain  terms  are  the  only  terms  that  are  not  zero 

at  the  wall  for  Uv,  with  values  larger  than  the  peak  values  observed  within  the  layer.  Both 
terms  are  reduced  to  near  zero  values  at  y^=5  before  they  reach  maxima  values  at  the 
wall.  They  are  approximately  equal  to  each  other  with  an  opposite  sign  near  the  wall 
below  y=5.  Away  from  the  wall  pressure-strain  and  production  terms  are  the  dominant 
terms  on  the  injection  side.  On  the  suction  side  the  pressure-diffusion  term  is  also  as  high 
as  these  two  terms.  Pressure-diffusion  changes  sign  three  times,  at  y^’^S,  near  y"=1 0  and 
one  around  y*«35.  Pressure-strain  is  positive.  Turbulent  diffusion  is  also  high  for  both 
cases.  These  trends  are  also  very  similar  to  the  2-D,  Station  5,  and  Separation  station 
data  of  the  present  study. 

V-B)  Study  of  Perot-Moin  (1993) 

Perot  studied  cases  of  Rej=—=54  and  134  in  a  temporally  developing  boundary 

ev 

layer  obtained  by  instantaneously  inserting  a  moving  solid  wall  to  an  isotropic 
homogenous  turbulence  (1993).  The  development  of  the  layer  depends  on  time  and 
distance  away  from  the  wall.  The  solutions  were  obtained  in  time  spanning  to  3  large 

eddy  turnover  time.  ( TQ=—=\avge  eddy  turnover  at  initial  conditions).  The  turbulence  is 

£ 

decaying  over  time.  The  flow  can  be  thought  as  a  2-D  TBL  with  high  free-stream 
turbulence. 

DNS  data  for  IP  show  that  the  production  and  diffusion  peak  at  closeby  distances 

away  from  the  wall  for  ReT=54  case  (Uo/L/^g„=0.5,  Uo=initial  rms  of  u  fluctuations). 
Pressure-diffusion  is  close  to  zero.  These  observations  are  similar  to  the  present  2-D 
data.  However,  the  following  observations  show  marked  differences  between  data  and 
DNS  solutions:  Dissipation  rate  is  the  highest  term.  Pressure-strain  is  as  high  as  the 
production  term.  Turbulent-diffusion  is  positive  near  the  wall  and  changes  sign  away  from 
the  wall. 

For  P  the  DNS  shows  that  pressure-diffusion  is  observed  as  high  as  the  pressure- 
strain  near  the  wall,  and  dissipation  rate  is  larger  than  both  of  the  terms  away  from  the 
wall.  Viscous-diffusion  increased  near  the  wall  and  reduced  to  zero  at  the  wall.  Production 

is  approximately  zero.  Similar  observations  were  made  for  the  2-D  flow  P  transport 
budget  terms  with  the  anisotropic  dissipation  rate  approximation.  DNS  solutions  show  that 
the  pressure-diffusion  term  is  one  of  the  most  important  terms  shaping  the  structure  of 
the  near-wall  flow.  Turbulent  transport  was  positive  throughout  the  layer,  which  is  not 
supported  by  the  present  data. 

At  the  wall  viscous-diffusion  and  dissipation  rate  terms  are  the  only  terms  different 
than  zero  for  vP.  Convection  and  pressure-diffusion  terms  are  practically  zero.  These  are 
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the  only  observations  from  the  2-D  flow  data  that  agree  with  the  DNS  solutions.  DNS 
shows  that  dissipation  rate  is  higher  than  pressure-strain  and  that  viscous  diffusion  is  as 
high  as  the  pressure-strain  near  the  wall. 

Pressure-strain  and  the  production  terms  are  the  largest  contributors  to  the 

transport  rate  of  Uv  with  opposite  signs,  which  are  similar  obsen/ations  from  the  present 
2-D  data.  The  following  observations  are  not  made  in  the  present  2-D  data.  Dissipation 
rate  and  viscous-diffusion  are  approximately  equal  to  each  other  with  opposite  sign.  The 
pressure-diffusion  term  is  zero.  Turbulent  diffusion  is  positive  near  the  wall  and  changes 
sign  to  show  a  negative  peak  at  the  location  where  production  profile  peaks. 

V-C)  Study  of  Mansour-Kim-Moin  (1988) 

Mansour,  Kim,  and  Moin  studied  Reynolds-stress  and  dissipation-rate  budgets  in 
a  turbulent  channel  flow  using  DNS.  The  Reynolds  number  of  the  flow  was 
Re^=^5/v=180  (5=channel  half  width).  Pressure-velocity  gradient  terms  were  not  split  into 
pressure-strain  and  pressure-diffusion,  but  together. 

For  IF,  variations  of  the  DNS  data  and  the  present  2-D  data  are  in  good 
agreement  globally.  For  the  low  Re  number  channel  flow  the  terms  are  closer  to  zero  at 
y'^ISO  where  in  the  present  data  the  values  are  different  than  zero  till  y"=1500.  Also  the 
scales  of  the  variations,  even  though  nondimensionalized  with  the  same  parameter,  are 
different  from  each  other. 

For  v^,  TP,  Uv,  agreement  between  the  data  and  the  DNS  simulations  are  globally 
very  good.  The  variations  of  the  dissipation  rate,  turbulent  diffusion  and  viscous-diffusion 
terms  within  the  layer  are  in  very  good  agreement  with  the  2-D  data. 

The  TKE  budget  term  variations  of  the  DNS  and  the  2-D  data  are  in  good 
agreement.  DNS  data  show  that  the  pressure-diffusion  term  is  very  close  to  zero,  except 
below  y'^^IO  with  a  zero  value  at  the  wall.  Viscous-diffusion  and  dissipation  rate  are  the 
only  terms  that  are  different  than  zero  at  the  wall  and  their  peak  values  occur  at  the  wall. 
The  differences  obsen/ed  are:  the  pressure-diffusion  term  of  the  present  data  is  slightly 
higher  than  the  DNS  data  relatively,  the  negative  peak  values  obtained  for  the  viscous- 
diffusion  and  the  turbulent-diffusion  are  not  at  the  same  level  for  the  2-D  data  as 
obsen/ed  in  DNS.  The  dip  of  the  dissipation  rate  profile  below  y"=15  is  not  observed  for 
the  DNS  data. 

They  analyzed  the  near  wall  asymptotic  behavior  of  the  Reynolds  stress  transport 

terms.  For  IP,  the  viscous-diffusion  term  becomes  positive  nearer  to  the  wall  than  the 
turbulent-diffusion  term.  Viscous-diffusion  at  the  wall  is  a  positive  constant.  The  pressure- 
gradient  term  asymptotically  varies  with  p  near  the  wall  and  is  zero  at  the  wall. 
Dissipation  rate  is  a  positive  constant  at  the  wall. 

For  P,  all  dissipation  rate,  viscous-diffusion,  and  pressure-gradient  terms  vary  as 
towards  the  wall.  Variations  dictate  that  while  the  diffusion  and  dissipation  rates  near 
the  wall  should  vary  as  positive,  pressure-gradient  term  close  to  the  wall  should  be 
negative.  This  is  observed  to  be  true  for  the  present  case  only  with  the  anisotropic 
dissipation  rate  approximation. 

Viscous-diffusion  and  dissipation  rate  terms  for  W  transport  are  positive  constants 
at  the  wall.  The  pressure-gradient  term  should  be  positive  near  the  wall,  which  is  also 


observed  for  the  present  data  as  the  sum  of  the  pressure-strain  and  pressure-diffusion 

terms.  For  uv,  analysis  showed  that  while  dissipation  rate,  viscous-diffusion,  and 
pressure-gradient  vary  with  y*,  production  and  turbulent  diffusion  vary  with  y^. 

Their  further  analysis  by  splitting  the  pressure-gradient  term  as  sum  of  pressure- 
strain  and  pressure-diffusion  showed  that  the  pressure-strain  and  the  pressure-diffusion 
terms  show  peaks  very  close  to  the  wall  around  y*=5  and  are  both  zero  at  the  wall. 
Similar  variations  are  observed  for  the  2-D  data  with  an  anisotropic  dissipation  rate 
approximation  further  away  from  the  wall  in  the  y^=10  to  20  region.  This  also  confirms  the 
fact  that  the  pressure-diffusion  term  affects  the  near  wall  distribution  highly. 

V-D)  Study  of  Moln-Shih-Driver-Mansour  (1990) 

Moin,  Shih,  Driver,  and  Mansour  computed  a  three-dimensional  turbulent  boundary 
layer  (3DTBL)  generated  with  a  sudden  pressure  gradient  (9P/9z=10  x^/6)  on  a  2-d  fully 
developed  plane  channel  flow.  The  Reynolds  number  based  on  channel  half  width  was 
3300.  The  flow  is  a  temporally  evolving  flow  with  zero  streamwise  velocity  gradients  rather 
than  a  spatially  developing  pressure-driven  3DTBL. 

In  a  3DTBL  turbulent  kinetic  energy  (TKE)  was  reduced.  The  W  increased  whileTF 

and  decreased.  Townsend’s  A,  structural  parameter  decreased  with  three- 
dimensionality.  The  shear-stress  angle  lagged  behind  the  flow-gradient  angle  near  the 
wall  and  led  away  from  the  wall.  The  pressure-strain  was  reduced  in  general  for  the  3-D 
case. 

Increased  three-dimensionality  reduced  the  TKE.  Turbulent  transport  and  pressure- 
diffusion  are  reduced.  Dissipation  rate  is  increased.  The  production  term  was  observed 
as  the  first  term  to  respond  to  the  three-dimensionality. 

For  v^,  appreciable  reduction  in  the  velocity-pressure  gradient  term  was  observed. 

Production  of  Hv  decrease,  which  is  believed  to  be  due  to  the  reduction.  Energy  is 

transferred  from  TF  into  which  is  believed  to  be  due  to  the  flow  turning. 

A  chain  of  events  causing  the  decrease  of  TKE  was  described.  Three- 

dimensionality  reduces  the  pressure-strain  in  the  "F  equation,  which  in  turn  reduces  the 

values.  Reduction  of  reduces  the  Uv,  and  reduction  of  Uv  reduces  TKE  production, 
which  in  turn  reduces  the  TKE. 

In  comparison  the  present  study  shows  that  the  nondimensional  terms  of  the  TKE 
transport  at  the  2-D  station  and  at  Station  5  are  not  much  different  from  each  other  with 
slightly  elevated  magnitudes  at  Station  5,  especially  of  the  dissipation  rate  term.  However 
the  terms  decay  faster  away  from  the  wall  at  Station  5  as  compared  to  the  2-D  station. 

For  the  present  data  there  is  a  reduction  in  the  pressure-strain  term  for  the  'F  transport 
and  the  reduction  of  the  production  term  magnitude  for  the  uv  transport  at  Station  5  as 

compared  to  the  2-D  station.  The  IF  and  Uv  data  also  show  the  reduction  of  these 
quantities  with  three-dimensionality  at  Station  5  in  nondimensional  form. 

V-E)  Study  of  Spalart  (1989) 
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Spalart  studied  a  3-D  boundary  layer  on  an  infinite  plate  generated  by  the  steadily 
rotating  free-stream  velocity  for  three  different  Reynolds  numbers.  The  Reynolds  numbers 
of  the  flow  based  on  the  laminar  boundary  layer  thickness  studied  are 

Re,=L/o5,/v=500,620,  and  767  If ,  f=angular  rate,  Uo=free-stream  velocity). 

Rotation  had  little  effect  on  the  turbulence.  The  viscous-diffusion  plus  dissipation  rate  and 
convection  plus  production  terms  along  with  the  pressure-gradient  term  and  a  time- 
dependent  term  were  presented.  Turbulent-diffusion  is  not  presented. 

The  TKE  and  Reynolds  stress  budgets  are  presented  for  Re|=767  case.  The  TKE 
distribution  shows  that  production  is  equal  to  dissipation  rate  with  peaks  located  around 
y'^^S.  Pressure-diffusion  is  an  important  contributor  near  the  wall  below  y'"«15  and 
changes  sign  to  positive  away  from  the  wall.  In  the  present  study  both  at  Station  5  and 
the  2-D  station,  dissipation  rate  and  production  values  peak  at  y'^^IS  with  different  values 
of  the  peak  and  pressure-diffusion  term  is  important  up  to  y''«40.  As  pointed  out  by 
Spalart  due  to  the  low  Reynolds  number  of  the  study  the  flow  simulated  is  weakly  three- 
dimensional. 

The  IP  budget  shows  similar  trends  of  the  TKE  distribution  with  highly  elevated 
values  of  the  combined  pressure-strain  and  pressure-diffusion  term  and  lower  values  of 
the  convection  plus  production.  However  below  y‘"=50  viscous-diffusion  plus  dissipation 
rate  and  production  plus  convection  terms  are  larger  in  magnitude  than  the  pressure- 
strain  plus  pressure-diffusion  term.  In  the  present  study  for  the  2-D  flow  production  plus 
convection  is  3  times  larger  than  the  dissipation  rate  and  pressure-strain  plus  pressure- 
diffusion  term.  At  Station  5  the  dissipation  rate  term  is  closer  to  the  production  than  the 
2-D  station.  Present  data  show  similar  trends  to  each  other  with  or  without  the  anisotropic 
dissipation  rate  approximation. 

From  the  DNS  the  TP  budget  shows  that  the  pressure  term  is  positive  and  higher 
in  magnitude  than  the  dissipation  rate  plus  viscous-diffusion.  The  convection  plus 
production  term  changes  sign  from  positive  to  negative  around  y^~20.  Dissipation  rate 
plus  viscous-diffusion  term  changes  sign  to  positive  around  y*»5.  Similar  observations  are 
made  for  the  present  data  with  anisotropic  dissipation  rate  approximation.  For  the  2-D 
station,  data  show  that  convection  does  not  change  sign  but  is  close  to  zero  at  y^^^lOO; 
viscous  terms  change  sign  below  y=10.  At  Station  5  convection  plus  production  changes 
sign  around  y'"=40  and,  viscous  terms  become  positive  around  y"=5.  Pressure  terms  are 
close  to  but  larger  than  the  dissipation  rate  plus  viscous-diffusion  in  magnitude. 

The  DNS  budget  for  TiP  shows  that  pressure  terms  are  lower  in  magnitude  than 
the  viscosity  involving  terms.  Production  plus  convection  is  negative  near  the  wall  below 
y'^IS.  Pressure  terms  change  sign  from  positive  to  negative  away  from  the  wall  at 
y/5=0.1.  At  Station  5  with  or  without  the  anisotropic  dissipation  rate  approximation,  data 
show  that  pressure  terms  are  low  in  magnitude  than  the  viscous  terms;  however,  the 
pressure  term  is  negative  below  y^^lS  opposite  to  the  trend  in  the  DNS,  it  becomes 
positive  between  y'^^IS  to  40  and  negative  above  this  y*  value.  The  production  plus 
convection  term  does  not  change  sign  as  calculated  by  DNS. 

The  sum  of  the  pressure  terms  in  the  -Uv  DNS  budget  are  close  in  magnitude  to 
the  production  plus  convection  terms.  The  dissipation  rate  term  changes  sign  to  negative 
at  y*=5.  Away  from  the  wall  above  y/5>0.35  production  plus  convection  term  and  pressure 
terms  change  sign.  The  dissipation  rate  term  is  comparatively  very  small.  For  the 
presently  reported  data  pressure  terms  and  production  plus  convection  terms  follow  each 
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other  in  magnitude  closely  similar  to  DNS  results,  however  they  do  not  change  sign  away 
from  the  wall. 

The  -  viv  budget  shows  that  the  pressure  terms  and  the  production  plus  convection 
terms  follow  each  other  closely  in  magnitude.  Dissipation  rate  plus  viscous-diffusion  is 
relatively  small  compared  to  these  terms  above  y^=10.  Similarly,  pressure  terms  and  the 
production  term  at  Station  5  are  close  in  magnitude  to  one  another.  Dissipation  rate  and 
convection  terms  are  negligibly  small  throughout  the  layer.  Pressure  terms  and  production 
terms  change  sign  above  y^«100  for  the  present  data  which  is  not  observed  for  the  DNS. 

The  -Tiw  budget  shows  that  viscous  terms  show  a  negative  peak  at  y'"«10  and 
pressure  terms  around  y*=25.  Positive  peak  of  production+convection  is  at  y'^IS.  All  the 
peak  values  are  comparable  to  each  other.  Dissipation  rate  decay  is  faster  than  the  other 
terms.  Present  data  with  anisotropic  dissipation  rate  approximation  show  that  all  terms 
in  comparison  to  DNS  peak  at  the  same  y^  location  with  higher  peak  magnitude  value  of 
production.  The  data  show  that  pressure  term  and  production  term  change  sign  away 
from  the  wall  but  at  a  lower  y*  location  then  the  DNS  results  show. 
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CHAPTER  VI  -  CONCLUSIONS 


Some  conclusions  of  the  present  study  can  be  made  about  the  effect  of  certain 
terms  in  the  transport  equations: 

1)  The  pressure-diffusion  term  largely  affects  the  transport  rate  budget  of  thelT^ 

normal  stress,  and  Uv  and  vw  shear  stresses  at  all  stations.  The  effect  of  the  term  is 
negligible  for  the  transport  rate  budget  of  the  TKE  and  other  stresses. 

2)  At  Station  5,  Separation  and  Vortex-core  stations  the  distributions  of  the 
transport  budget  terms  decrease  progressively  more  steeply  to  approximately  zero  values 

than  the  2-D  flow  station  at  lower  y"^  locations.  This  is  more  clearly  visible  for  the  IF,  IF, 

and  Uv  stresses.  At  the  Vortex  station  these  observations  are  accompanied  with  the 
effects  of  the  outer  region  vortex.  Vortex  affects  the  terms  to  be  away  from  zero.  At  the 

2-D  station  Uw  and  Uw  stresses  are  zero. 

3)  The  effects  of  not  using  a  near  wall  correction  and  an  isotropic  dissipation  rate 
approximation  were  examined.  Model  predictions  and  the  data  agreement  deteriorates 
with  the  increased  complexity  of  the  flow.  Very  near  wall  data  values  and  predictions  are 
different  from  each  other  for  all  stations  below  y*=35.  Overall  different  models  predict 
different  stress  transport  pressure-strain  terms  better.  The  FLT2  model  behaves  better 

in  predicting  the  IF,  and  stresses  for  complex  Station  5  and  Separation  stations  but 

not  the  other  stations.  The  TF  and  'F  estimates  are  very  different  than  the  data  for  the 

Vortex-core  station.  At  the  2-D  station  for  IF,  above  y'^=40  all  the  models  except  the  FLT 

models,  and  in  the  same  region  for  the  'F  all  the  models  except  the  FLT  and  OP  models 

work  best.  At  the  2-D  station  near  wall  variation  of  UF  is  well  approximated  by  FLT1 
model.  Above  y'^^SS,  SSG,  SLCL,  GL,  and  LRR  models  work  best.  At  Station  5,  FLT2 

model  is  better  above  y'^=30.  The  vP  pressure-strain  term  is  not  well  predicted  for 
Separation  and  Vortex-core  stations  by  any  model.  None  of  the  models  predicts  the 

pressure-strain  variation  of  any  shear  stress  transport  for  the  Vortex-core  station.  Them? 
stress  is  better  predicted  by  the  SLCL  model  for  2-D  and  Separation  stations.  For  Station 

5  LRR  and  SSG  models  predict  results  closer  to  the  data.  The  Uw  pressure-strain  value 

is  better  predicted  by  the  FLT1  model  and  vw  is  better  predicted  by  the  SLCL  models. 

4)  Anisotropic  dissipation  rate  effects  were  also  examined.  The  isotropic  dissipation 
rate  assumption  results  in  equally  distributed  dissipation  rate  for  the  normal  stresses  and 
no  dissipation  rate  for  the  shear  stresses.  On  the  other  hand  the  anisotropic  dissipation 
rate  model  of  Hallback  et.  al.  shows  that  the  dissipation  rate  is  increasingly  anisotropic 

with  increased  three-dimensionality  effects.  The  dissipation  rate  in  the  'F  normal  stress 

is  relatively  smaller  then  the  TF  and  UF  dissipation  rate  terms  for  the  3-D  flow  stations 

than  the  2-D  fiow  station.  Dissipation  rate  in  Uv  and  vw  stresses  are  on  the  order  of  half 

of  the  'F  terms.  The  Uw  dissipation  rate  near  the  wall  is  larger  than  the  'F  dissipation 
rate  with  an  opposite  sign  for  the  Station  5,  and  Separation  stations.  For  3-D  flow  stations 
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dissipation  rate  distributions  are  approximately  isotropic  above  y*=100.  However  for  2-D 
flow  it  becomes  isotropic  above  y*=350. 

Anisotropy  of  the  dissipation  rate  highly  changes  the  experimental  pressure-strain 

distributions  extracted  for  the  normal  stresses;  especially  of  the  normal  stress.  Near 

the  wall  between  y"  of  10-15  to  30-40,  pressure-strain  distributions  become  negative 
and  show  negative  peaks.  The  effect  of  an  anisotropic  distribution  for  the  shear  stress 
pressure-strain  transport  terms  is  non-negligible  but  small. 

Anisotropic  dissipation  rate  approximation  increases  the  differences  between  the 
predictions  and  the  data  especially  near  the  wall  below  y'^SO  for  all  the  models  except 
the  FLT2  model.  With  the  anisotropic  dissipation  rate  approximation  FLT2  model  works 

well  near  the  wall  especially  for  the  "F  for  2-D,  Station  5,  and  Separation  stations.  The 
models  which  work  best  for  the  other  stresses  with  isotropic  dissipation  rate 
approximation  also  work  best  with  the  anisotropic  dissipation  rate  approximation. 

Comparison  of  the  present  data  to  the  DNS  results  shows  that  the  present  data 
are  in  better  agreement  with  the  DNS  solutions  with  the  anisotropic  dissipation  rate 
approximation. 

5)  The  near-wall  corrections  with  the  isotropic  dissipation  rate  approximation  bring 
the  data  and  the  estimates  close  to  each  other  for  the  normal  stresses.  Corrections  for 

the  Uv  stresses  very  near  the  wall  are  very  large  and  corrections  to  the  Uw  data  do  not 

make  the  agreement  better,  vw  corrections  very  near  the  wall  below  y'^^IS  is  opposite 
to  the  trend  of  the  data.  Near  wall  corrections  result  in  LRR,  and  GL  models  to  be  much 
closer  to  the  data  especially  for  the  shear  stress  pressure-strain  values,  above  y'^IS. 

6)  The  combined  effect  of  near-wall  correction  and  the  anisotropic  dissipation  rate 

approximation  in  general  reduces  the  agreement  between  the  data  and  the  estimates  for 77^ 

and  v^.  The  agreement  is  not  affected  by  the  approximations.  Agreement  between 
the  shear  stress  data  and  predictions  with  the  combined  effect  is  similar  to  the  case  with 
only  the  near  wall  correction. 

V-A)  Future  work 

The  TKE  and  shear  stress  transport  budget  balance  were  extracted  using  the 
pressure-diffusion  model  of  Lumley  (1978).  Experimental  dissipation  rate  terms  extracted 
for  the  TKE  transport  equation  are  not  highly  affected  by  this  modelling.  However,  the 
pressure-strain  terms  extracted  from  the  shear  stress  transport  equations  show  that 
pressure-diffusion  is  one  of  the  most  prominent  terms  in  the  balance  using  Lumley’s 

approach,  especially  for  the  v^,  Uv,  and  vw  stresses.  The  validity  of  this  model  or  more 
improved  models  are  needed.  A  recent  study  by  Kim  and  Chung  (1995)  was  directed 
towards  this  issue. 

Other  important  needed  research  is  the  further  validation  of  the  anisotropic 
dissipation  rate  models.  The  dissipation  rate  term  (e),  extracted  from  the  TKE  budget  is 
distributed  among  the  stresses.  An  isotropic  dissipation  rate  approximation  results  in 
equally  distributed  dissipation  rate  terms  among  the  normal  stresses  and  no  dissipation 
rate  term  in  the  shear  stresses.  Anisotropic  dissipation  rate  models  which  are 
asymptotically  correct  near  the  wall  did  not  seem  to  vary  differently  then  the  isotropic 
distributions  for  the  high  Ree  cases  studied  in  the  present  report  away  from  the  wall.  The 
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Hallback,  Groth,  and  Johansson  anisotropic  dissipation  rate  model  not  only  satisfied  near 
wall  variations  but  also  resulted  in  different  variations  throughout  the  layers  for  different 
stresses.  Both  isotropic  dissipation  rate  and  anisotropic  dissipation  rate  approximations 
were  studied  in  this  report.  Presented  budget  distributions  with  anisotropic  dissipation  rate 
approximations  are  closer  to  the  DNS  results,  which  is  a  further  support  for  the 
anisotropic  dissipation  rate  model  chosen  in  the  present  study.  However,  the  Reynolds 
numbers  of  the  DNS  are  usually  very  low  to  make  complete  direct  comparisons. 
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Tables 


TABLE.2  SOME  LENGTH  SCALES  OBTAINED  FROM  LDV  DATA  IN  TUNNEL 
COORDINATES. 

Stations 

8 

(mm) 

5, 

(mm) 

5^ 

(mm) 

53 

(mm) 

84 

(mm) 

85 

(mm) 

83 

(mm) 

87 

(mm) 

2-D 

station 

39.1 

6.28 

6.27 

4.55 

-0.29 

0.04 

-0.26 

-0.00 

Station  5 

40.0 

6.26 

5.38 

4.64 

7.47 

-1.29 

6.18 

-1.39 

Separatio 
n  station 

39.3 

6.40 

5.16 

4.71 

9.04 

-1.54 

7.49 

-1.96 

Vortex 

station 

40.9 

6.59 

4.78 

4.85 

10.67 

-1.92 

8.75 

-2.83 

Ug=Velocity  magnitude  at  the  layer  edge. 


5=  y  where  yCF+W^ I Ug=0. 995=  Boundary  Layer  Thickness 

OO 

5^  =  [(1 --^)dy=  Streamwise  Displacement  Thickness 

'i  u 

0  e 

5a=  [(1  )dy=  Magnitude  Displacement  Thickness 

0 

OO 

83=  [(1 --^)-^c/y=  Streamwise  Momentum  Thickness 

u  u 

0  e 

64=  r(-J^)afy=  Lateral  Displacement  Thickness 

«  u 

0  e 

^^=j■^{^-J:L)dy=  Lateral  Momentum  Thickness 

5g=J-Jl^!^dy=  Cross  Product  Momentum  Thickness 

0 

7 

57=1- — -dy=  Crossflow  Momentum  Thickness 
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Figure  2.  Oil-flow  visualization  picture  obtained  on  the  tunnel  floor. 
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Figure  3.  Velocity  magnitude  contours  and  secondary  flow  streamlines  measured  at  stations  along  the  23°  line. 
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Figure  6.  W/u^  mean  velocity  at  different  stations  in  tunnel  coordinates.  •,  2-D;  O,  Station  5,  ■,  Separation  station; 
Vortex  station. 
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Figure  7.  IFlu^  normal  stress  component  at  different  stations  in  tunnel  coordinates.  2-D;  O,  Station  5,  ■,  Separation 
station;  □,  Vortex  station. 
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Figure  8.  normal  stress  component  at  different  stations  in  tunnel  coordinates.  •,  2-D;  O,  Station  5,  ■,  Separation 
Station;  □,  Vortex  station. 


Figure  9.  w'^/u^  normal  stress  component  at  different  stations  in  tunnel  coordinates.  2’D,  O,  Station  5,  H,  Separation 
station;  □,  Vortex  station. 
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Figure  1 0.  uv/uf  shear  stress  component  at  different  stations  in  tunnel  coordinates.  •,  2-D:  O,  Station  5,  ■,  Separation 
station;  □,  Vortex  station. 


Figure  1 1 . 1ml shear  stress  component  at  different  stations  in  tunnel  coordinates.  •,  2-D;  O,  Station  5,  ■,  Separation 
station;  □,  Vortex  station. 
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Figure  14.  TFwIu^  triple  product  at  different  stations  in  tunnel  coordinates.  #,  2-D;  O,  Station  5,  ■,  Separation  station; 
□,  Vortex  station. 


Figure  15.  v^w/u^  triple  product  at  different  stations  In  tunnel  coordinates.  •,  2-0;  O,  Station  5,  ■,  Separation  station; 
□,  Vortex  station. 
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Figure  16.  IFulu^  triple  product  at  different  stations  in  tunnel  coordinates.  •,  2-D:  O,  Station  5,  ■>  Separation  station; 
□,  Vortex  station. 


Figure  17.  Wulu^  triple  product  at  different  stations  in  tunnel  coordinates.  •,  2-D;  O,  Station  5,  ■,  Separation  station; 
□,  Vortex  station. 
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Figure  18.  "Wvlu^  triple  product  at  different  stations  in  tunnel  coordinates.  •,  2-D;  O,  Station  5,  ■,  Separation  station; 
□,  Vortex  station. 


Figure  19.  Uvwlu^  triple  product  at  different  stations  in  tunnel  coordinates.  •,  2-D;  O,  Station  5,  ■,  Separation  station; 
□,  Vortex  station. 
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Figure  20.  TFlu^  triple  product  at  different  stations  in  tunnel  coordinates.  •,  2-D;  O,  Station  5,  ■,  Separation  station,  □, 
Vortex  station. 


Figure  21 .  v^/u^  triple  product  at  different  stations  In  tunnel  coordinates.  •,  2-D;  O,  Station  5,  ■.  Separation  station;  □, 
Vortex  station. 
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Figure  24.  Flow  gradient  angle  (FGA)  calculated  at  different  stations  in  tunnel  coordinates.  •,  2-0;  O,  Station  5, 
Separation  station;  □,  Vortex  station. 
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Figure  25.  Shear  stress  angle  (SSA)  calculated  at  different  stations  in  tunnel  coordinates.  •,  2-D;  O,  Station  5, 
Separation  station;  □,  Vortex  station. 


Figure  29.  TF  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  2-D  station  in  tunnel  coordinates. 
•,  convection  (C);  O,  production  (PR);  0,  pressure-strain  ((}));  □,  turbulent  diffusion  (TD);  pressure  diffusion  (PD);  v, 
viscous  diffusion  (VD);  ♦,  dissipation  (e). 


Figure  30.  "F  normal  stress  transport  budget  with  Isotropic  dissipation  assumption  at  2-D  station  in  tunnel  coordinates. 
•,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (([));  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion  (PD);  v, 
viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  31 .  vP  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  2-D  station  in  tunnel  coordinates. 
•,  convection  (C);  O,  production  (PR);  0,  pressure-strain  ((());  □,  turbulent  diffusion  (TD);  *■,  pressure  diffusion  (PD);  v, 
viscous  diffusion  (VD);  ♦,  dissipation  (e). 


Figure  32.  W  shear  stress  transport  budget  at  2-D  station  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR); 
0,  pressure-strain  (<!>);  □,  turbulent  diffusion  (TD);  pressure  diffusion  (PD);  v.  viscous  diffusion  (VD).  Dissipation  is 
assumed  zero. 
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Figure  33.  uw  shear  stress  transport  budget  at  2-D  station  in  tunnel  coordinates.  #,  convection  (C);  O,  production  (PR); 
0,  pressure-strain  ((|));  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation  is 
assumed  zero. 
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Figure  34.  vw  shear  stress  transport  budget  at  2-D  station  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR); 
0,  pressure-strain  ((j));  □,  turbulent  diffusion  (TD);  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation  is 
assumed  zero. 
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Figure  35.  Turbulent  kinetic  energy  transport  budget  at  2-D  station  in  tunnel  coordinates.  •,  convection  (C);  O, 
production  (PR);  □,  turbulent  diffusion  (TD);  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
Pressure-strain  term  is  zero  by  definition. 
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Figure  36.  IF  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Station  5  in  tunnel  coordinates. 
•,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (^);  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion  (PD);  v, 
viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  37.  "F  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Station  5  in  tunnel  coordinates. 
#,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (([));  □,  turbulent  diffusion  (TD);  pressure  diffusion  (PD);  v, 
viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  38.  IF  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Station  5  in  tunnel  coordinates. 
•,  convection  (C);  O,  production  (PR);  0,  pressure-strain  □,  turbulent  diffusion  (TD);  *•,  pressure  diffusion  (PD);  v, 
viscous  diffusion  (VD);  ♦,  dissipation  (e). 


Figure  39.  Uv  shear  stress  transport  budget  at  Station  5  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR), 
0,  pressure-strain  {^Y  □.  turbulent  diffusion  (TD);  *.  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation  is 
assumed  zero. 
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Figure  40.  Uw  shear  stress  transport  budget  at  Station  5  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR); 
0,  pressure-strain  (cj)):  □,  turbulent  diffusion  (TD);  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation  is 
assumed  zero. 
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Figure  41.  vw  shear  stress  transport  budget  at  Station  5  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR); 
0,  pressure-strain  (([));  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation  is 
assumed  zero. 
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Figure  42.  Turbulent  kinetic  energy  transport  budget  at  Station  5  in  tunnel  coordinates.  •,  convection  (C);  O,  production 
(PR):  □,  turbulent  diffusion  (TD);  *,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e).  Pressure-strain 
term  is  zero  by  definition. 


Figure  43.  TF  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Separation  station  in  tunnel 
coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  ((|));  □,  turbulent  diffusion  (TD);  pressure  diffusion 
(PD);  V,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  44.  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Separation  station  in  tunnel 
coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (([));  □,  turbulent  diffusion  (TD);  pressure  diffusion 
(PD);  V,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  45.  W  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Separation  station  in  tunnel 
coordinates.  •,  convection  (C);  O.  production  (PR);  0.  pressure-strain  (<(>);  □,  turbulent  dittusion  (TD);  pressure  diffusion 
(PD):  V,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 


Figure  46.  Uv  shear  stress  transport  budget  at  Separation  station  in  tunnel  coordinates.  •,  convection  (C);  O,  production 
(PR);  0.  pressure-strain  (4.):  □.  turbulent  diffusion  (TD);  pressure  diffusion  (PD);  v.  viscous  diffusion  (VD).  Dissipation 
Is  assumed  zero. 
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Figure  47.  uw  shear  stress  transport  budget  at  Separation  station  in  tunnel  coordinates.  •,  convection  (C);  O,  production 
(PR);  0,  pressure-strain  ((}));  □,  turbulent  diffusion  (TD);  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation 


is  assumed  zero. 
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Figure  48.  vw  shear  stress  transport  budget  at  Separation  station  in  tunnel  coordinates.  •,  convection  (C);  O,  production 
(PR);  0,  pressure-strain  (<)));  turbulent  diffusion  (TD);  a,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation 
is  assumed  zero. 
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Figure  49.  Turbulent  kinetic  energy  transport  budget  at  Separation  station  in  tunnel  coordinates.  •,  convection  (C);  O, 
production  (PR);  □,  turbulent  diffusion  (TD);  a.  pressure  diffusion  (PD);  v.  viscous  diffusion  (VD);  ♦,  dissipation  (e) 
Pressure-strain  term  is  zero  by  definition. 


Figure  50.  7F  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Vortex  station  In  tunnel 
coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (<t));  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion 
(PD);  V,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  51.  'F  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Vortex  station  in  tunnel 
coordinates.  #,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (<[));  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion 
(PD);  V,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  52.  TF  normal  stress  transport  budget  with  isotropic  dissipation  assumption  at  Vortex  station  in  tunnel 
coordinates.  •.  convection  (C);  O,  production  (PR);  0,  pressure-strain  W;  □,  turbulent  diffusion  (TD);  pressure  diffusion 
(PD);  V,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 


Figure  53.  uv  shear  stress  transport  budget  at  Vortex  station  in  tunnel  coordinates.  •,  convection  (C);  O,  production 
(PR);  0,  pressure-strain  (()>);  □,  turbulent  diffusion  (TD);  *,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation 
is  assumed  zero. 
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Figure  54.  Uw  shear  stress  transport  budget  at  Vortex  station  in  tunnel  coordinates.  #,  convection  (C);  O,  production 
(PR);  0,  pressure-strain  (([));  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation 
is  assumed  zero. 


Figure  55.  vw  shear  stress  transport  budget  at  Vortex  station  in  tunnel  coordinates.  •,  convection  (C);  O,  production 
(PR);  0,  pressure-strain  (([));  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD).  Dissipation 
is  assumed  zero. 
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Figure  56.  Turbulent  kinetic  energy  transport  budget  at  Vortex  station  in  tunnel  coordinates.  •,  convection  (C);  O, 
production  (PR);  □,  turbulent  diffusion  (TD);  a,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
Pressure-strain  term  is  zero  by  definition. 
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Figure  57.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  »;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  58.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  59.  Comparison  of  W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 


-0-* 
-•-LRR  SB 
-e-LRR  SL 
-■-OP 
-A-SSG 
-^LCL 
-iSr-fLTSB 
-'F-FLTSL 

UV,  TO 
2-D 


Figure  60.  Comparison  of  UP  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 


Figures 


84 


y" 


Figure  61 .  Comparison  of  'uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  62.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  63.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 


Figure  64.  Comparison  of  "F  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale-Sarkar-Gatski 
(SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  65.  Comparison  of  W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  66.  Comparison  of  uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  67.  Comparison  of  Uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strairi 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  68.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  •  Speziale-Sarkar-Gatski 
(SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  69.  Comparison  of  W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  70.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale- 
Sarkar-Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  71 .  Comparison  of  TP  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  72.  Comparison  of  Tiv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  73.  Comparison  of  Uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  74.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  75.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  76.  Comparison  of  "F  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  77.  Comparison  of  "W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  Hi  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  isotropic. 
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Figure  78.  Comparison  of  Vv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  79.  Comparison  of  Uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 
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Figure  80.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  dissipation  is  assumed  zero. 


Figures 


94 


Anisotropic  dissipation  tensor  of  Hallback,  Groth,  and  Johansson,  without  the  near  wall  corrections  of  Launder-Reece- 
Rodi.  The  balance  figures  change  since  the  dissipation  and  pressure-strain  terms  are  effected  from  these  anisotropies. 
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Figure  81.  Anisotropic  dissipation  rate  distribution  calculated  using  Hallback  et  al.  (1990)  model  for  the  2-D  station  in 
tunnel  coordinates.  IF,  •;  "F,  O;  FF,  ■;  W,  □;  TTiv,  w),  isotropic  distribution,  +. 
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Figure  82.  Anisotropic  dissipation  rate  distribution  calculated  using  Hallback  et  al.  (1990)  model  for  Station  5  in  tunnel 
coordinates.  TF,  •;  F,  O;  W,  ■;  Fv,  □;  uw,  Fw,  isotropic  distribution,  +. 
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Figure  83.  Anisotropic  dissipation  rate  distribution  calculated  using  Hallback  et  al.  (1990)  model  for  Separation  station 
in  tunnel  coordinates.  IF,  •;  'F,  O;  vF,  ■;  uv,  □;  "Uw,  Vw,  ▼;  isotropic  distribution,  +. 
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Figure  84.  Anisotropic  dissipation  rate  distribution  calculated  using  Hallback  et  al.  (1990)  model  for  Vortex  station  in 
tunnel  coordinates.  TF,  •;  "F,  O;  W,  ■;  uv,  □;  Fw,  vw,  ▼;  isotropic  distribution,  +. 
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Figure  85.  TF  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1 990)  at  2-D  station 
in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (([)):  □,  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 


Figure  86.  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  2-D  station 
in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (<|));  □,  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  87.  W  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  2-D  station 
in  tunnel  coordinates.  •,  convection  (C);  O.  production  (PR);  0,  pressure-strain  («l>);  □.  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  88.  TTv  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  2-D  station 
in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0.  pressure-strain  (().);  □,  turbulent  diffusion  (TD);  a, 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  89.  IP  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Station  5 
in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (((•);  □,  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  90.  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Station  5 
in  tunnel  coordinates.  #,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (({));  □,  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  91.  TF  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  HallbSck  et  al  (1990)  at  Station  5 
in  tunnel  coordinates.  •,  convection  (C);  O.  production  (PR):  0.  pressure-strain  (<))):  □,  turbulent  diffusion  (TD);  a, 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  92.  uv  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  HallbSck  et  al  (1990)  at  Station  5  in 
tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0.  pressure-strain  m  □,  turbulent  diffusion  (TD);  a,  pressure 
diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  93.  Uw  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Station  5  in 
tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (<[));  □,  turbulent  diffusion  (TD);  pressure 
diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  94.  Vw  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Station  5  in 
tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (^)\  □,  turbulent  diffusion  (TD);  a,  pressure 
diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  95.  IP  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Separation 
station  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (<|>):  □,  turbulent  diffusion  (TD); 
A,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  96.  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Separation 
station  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (<|));  □,  turbulent  diffusion  (TD); 
A,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  97.  TF  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Separation 
station  In  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (()));  □,  turbulent  diffusion  (TD); 
A,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  98.  Uv  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Separation 
station  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  ((|>);  □,  turbulent  diffusion  (TD); 
A,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  99.  Hw  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Separation 
station  in  tunnel  coordinates.  •.  convection  (C);  O,  production  (PR);  0,  pressure-strain  (())):  □,  turbulent  diffusion  (TD); 
A,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  100.  vw  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Separation 
station  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (^);  □,  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  101.  IF  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Vortex 
station  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  ((]));  □,  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  102.  IF  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Vortex 
station  in  tunnel  coordinates.  #,  convection  (C);  O,  production  (PR);  0,  pressure-strain  ((());  □,  turbulent  diffusion  (TD); 
A,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  103  vF  normal  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Vortex 
station  in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR):  0,  pressure-strain  (())):  □.  turbulent  diffusion  (TD); 
*,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  1 04.  TJv  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Vortex  station 
in  tunnel  coordinates.  •,  convection  (C);  O,  production  (PR);  0,  pressure-strain  ((]));  □.  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  105.  Vw  shear  stress  transport  budget  with  anisotropic  dissipation  mode!  of  Hallback  et  al  (1990)  at  Vortex 
station  in  tunnel  coordinates.  #,  convection  (C);  O,  production  (PR);  0,  pressure-strain  (<[));  □,  turbulent  diffusion  (TD); 
A,  pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 


Figure  106.  vw  shear  stress  transport  budget  with  anisotropic  dissipation  model  of  Hallback  et  al  (1990)  at  Vortex 
station  in  tunnel  coordinates.  #,  convection  (C);  O,  production  (PR);  0,  pressure-strain  ((!>);  □,  turbulent  diffusion  (TD); 
pressure  diffusion  (PD);  v,  viscous  diffusion  (VD);  ♦,  dissipation  (e). 
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Figure  107.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-RodI  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 08.  Comparison  of  "P  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 09.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  110.  Comparison  of  Uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  »;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  111.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 1 2.  Comparison  of  "F  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  m,  Speziale-Sarkar-Gatski 
(SSG),  *:  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wail  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 1 3.  Comparison  of  W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  114.  Comparison  of  W  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  115.  Comparison  of  uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2).  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 1 6.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 1 7.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih>Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 1 8.  Comparison  of  "P  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-RodI  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback, 
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Figure  119.  Comparison  of  vF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  No  near-wall  corrections,  anisotropic  dissipation  of  HallbSck. 
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Figure  120.  Comparison  of  Uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■,  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  121 .  Comparison  of  uiv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  122.  Comparison  of  viv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 23.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  124.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 25.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  126.  Comparison  of  uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  1 27.  Comparison  of  uiv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters(OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  {FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Figure  128.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters(OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  No  near-wall  corrections,  anisotropic  dissipation  of  Hallback. 
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Near  wall  correction  of  Launder-Reece-Rodi  for  their  model  as  in  Shih-Lumley  paper,  for  different  flow  conditions 
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Figure  129.  Near-wall  correction  to  the  pressure-strain  rate  distribution  calculated  using  Launder-Reece-Rodi  (1975) 
model  for  the  2-D  station  in  tunnel  coordinates.  IP,  •;  “P,  O;  W,  ^  Pv,  □;  Pw,  *;  vw,  ▼. 

pressure-strain  near  wall  correction  terms  for  3D-st5  by  LRRSL 
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Figure  130.  Near-wall  correction  to  the  pressure-strain  rate  distribution  calculated  using  Launder-Reece-Rodi  (1975) 
model  at  Station  5  in  tunnel  coordinates.  IP,  •;  P,  O:  TP,  ■;  Pv,  □:  Pw,  Tw,  ▼. 
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Figure  131.  Near-wall  correction  to  the  pressure-strain  rate  distribution  calculated  using  Launder-Reece-Rodi  (1975) 
model  for  Separation  station  in  tunnel  coordinates.  TF,  •;  "P,  O;  liF,  ■;  Uv,  □;  llw,  a;  w,  ▼. 

pressure— strain  near  wall  correction  terms  for  3D— vor  by  LRRSL 


Figure  132.  Near-wall  correction  to  the  pressure-strain  rate  distribution  calculated  using  Launder-Reece-Rodi  (1975) 
model  for  Vortex  station  in  tunnel  coordinates.  TF,  •;  IF,  O;  TiF,  ^  Uv,  □;  Uw,  a;  vw,  ▼. 
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Near  wall  correction  of  Launder-Reece-Rodi  Shih-Lumley  paper  applied  for  all  the  models 
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Figure  1 33.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 
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Figure  134.  Comparison  of  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 
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Figure  1 35.  Comparison  of  W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Gberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  K  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1 975),  dissipation  is  assumed  isotropic. 
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Figure  136.  Comparison  of  uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Figure  137,  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberiack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 


Figure  138.  Comparison  of  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberiack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 
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Figure  1 39.  Comparison  of  "W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumiey/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 
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Figure  140.  Comparison  of  w  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  *:  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Figure  141.  Comparison  of  Uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Figure  142.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Figure  1 43.  Comparison  of  IP  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  {FLT2),  v.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 
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Figure  1 44.  Comparison  of  "F  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 


I  I  I  ■TH — I  I  I  I  I  I  !"■  I  I  ■  I  ■  III  I  I  I  I  I  1 1  I  I  '  I  ■  I  M — 1  I  I  I  I  I  I  '  I  '  I’H'T 


# 


Figures 


126 


0.3 


h 

D 


E 


0) 


-•-LRRSB 

-e-LRRSL 

-■-OP 

-A-SSG 

-^-^LCL 

-A-FLTSB 

-9-FLTSL 

wMC 
3— D.sep 
near  wall 
correc. 
of  LRRSL 
applied 


y" 


Figure  1 45.  Comparison  of  'W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  K  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 
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Figure  1 46.  Comparison  of  Uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Figure  1 47.  Comparison  of  uiv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Figure  148.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 


Figures 


128 


0.20 

0.15 

0.10 

rA 

> 

h 

D 

0.05 

c 

c 

L. 

0) 

0.00 

-0.05 

-0.10 

«  »  ■  1  - 1  1  I  TTT 

»  I  M  * « ■  1 — 1  Mill 

■  '  •  1  ■  ‘  ■  1 

- 

i 

■ 

1 

«/(7/ 

\  1/ 

* 

- 

-#-URR  SB 

-e-iJRR  SL 

-■-OP 

-A-SSG 

-^SLCL 

-A-FLTSB 

-’^-FITSL 

u2,  TC 
3— D,vor 
near  wall 
correc. 
of  LRRSL 
applied 


100  2  3  101  2  3  102  2  3  1  03  2  3 


y" 


Figure  149.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack- Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 


SB 

-en^RR  SL 
-WHDP 

-▼-SLCL 

-^5r-FLTSB 

-^5M=LTSL 

vMC 
3— D,vor 
near  wall 
correc. 
of  LRRSL 
applied 


Figure  1 50.  Comparison  of  "P  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  m  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 
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Figure  1 51 .  Comparison  of  'wF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  isotropic. 
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Figure  1 52.  Comparison  of  uV  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■  Speziale-Sarkar- 
Gatskl  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Figure  153.  Comparison  of  uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters(OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Figure  154.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  dissipation  is  assumed  zero. 
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Both  near  wall  correction  of  LRR  and  anisotropic  dissipation  rate  of  Hallback  are  applied. 
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Figure  1 55.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL).  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 


Figure  1 56.  Comparison  of  "F  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  1 57.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station,  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  158.  Comparison  of  uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  the  2-D  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  1 59.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
{FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  1 60.  Comparison  of  "F  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR).  O;  Oberlack-Peters  (OP),  m  Speziale-Sarkar-Gatski 
(SSG),  *:  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1 975),  anisotropic  dissipation  of  Hallback  (1 990). 
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Figure  1 61 .  Comparison  of  W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  K  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tse!epidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  162.  Comparison  of  Uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5,  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  V  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  1 63.  Comparison  of  Uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar-Gatski 
(SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  164.  Comparison  of  w  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Station  5.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale-Sarkar-Gatski 
(SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis  2 
(FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  1 65.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder- Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback 
(1990). 


— 0— 

-■"LRR  SB 
-e-LRRSL 
-*-OP 
-A-SSG 

-▼-SLCL  ^ 

-A-FLTSB  • 

-V-FLTSL 

vMC 

3-D,sep 

anisotropi- 

dissipatior  # 

of 

Hallback 

and 

near— wall 
correct,  o 

LRRSL  • 


y" 


Figure  1 66.  Comparison  of  IF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback 
(1990). 
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Figure  1 67.  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  GIbson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  a;  Shih-Lumley/Chol-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback 
(1990). 
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Figure  1 68.  Comparison  of  Uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback 
(1990). 
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Figure  1 69.  Comparison  of  Hw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ^  Speziale- 
Sarkar-Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  v.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback 
(1990). 
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Figure  1 70.  Comparison  of  vw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Separation  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  m,  Speziale- 
Sarkar-Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder- 
Tselepidakis  2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback 
(1990). 
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Figure  171 .  Comparison  of  TF  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  1 72.  Comparison  of  'F  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  a;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  1 73.  Comparison  of  W  normal  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  174.  Comparison  of  uv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  175.  Comparison  of  'uw  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  *;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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Figure  176.  Comparison  of  viv  shear  stress  transport  experimental  pressure-strain  term  0,  to  theoretical  pressure-strain 
terms  for  Vortex  station.  Gibson-Launder,  •;  Launder-Reece-Rodi  (LRR),  O;  Oberlack-Peters  (OP),  ■;  Speziale-Sarkar- 
Gatski  (SSG),  *■;  Shih-Lumley/Choi-Lumley  (SLCL),  ▼;  Fu-Launder-Tselepidakis  1  (FLT1),  a;  Fu-Launder-Tselepidakis 
2  (FLT2),  V.  Near-wall  corrections  of  Launder-Reece-Rodi  (1975),  anisotropic  dissipation  of  Hallback  (1990). 
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